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Abstract 

This  paper  addresses  aspects  of  model  development  and  control  design  for  smart  structure  which 
utilize  piezoelectric,  electrostrictive,  magnetostrictive  or  shape  memory  alloys.  The  advent  of  these 
materials  as  sensors  and  actuators  provides  the  opportunity  for  designing  systems  with  improved  per¬ 
formance  while  reducing  weight,  hardware  and  power  requirements.  At  high  drive  levels,  however,  all 
of  these  materials  exhibit  nonlinear  constitutive  relations  and  varying  degrees  of  hysteresis  which  must 
be  accommodated  to  realize  their  full  potential  in  smart  structures.  We  outline  a  variety  of  modeling 
techniques  and  illustrate  their  use  in  feedback,  feedforward  and  adaptive  control  designs. 
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1  Introduction 


Increased  demands  for  high  performance  control  design  in  combination  with  recent  advances  in  material 
science  have  produced  a  class  of  systems  termed  smart,  intelligent  or  adaptive  structures.  While  subtle 
differences  may  be  associated  with  the  individual  terms,  these  structures  are  generally  defined  as  systems 
whose  dynamics  can  be  monitored  or  modified  by  distributed  sensors  and  actuators,  in  accordance  with 
an  integrated  control  law,  to  accommodate  time-varying  exogenous  inputs  or  changing  environmental 
conditions.  Specific  choices  for  the  actuators,  sensors  and  control  laws  are  dictated  by  the  design 
requirements  for  the  system. 

For  aeronautic  and  aerospace  systems,  control  transducers  must  be  light-  weight  and  should  typically 
have  minimal  effect  on  the  passive  system  dynamics.  Furthermore,  actuators  must  provide  the  required 
strain  or  force  inputs  using  the  available  power  supplies  which,  in  certain  aerospace  structures,  may 
require  the  scavenging  of  power  from  other  components  in  the  system.  Restrictions  on  size  and  weight 
also  dictate  that  transducers  in  some  regimes  must  be  capable  of  multiple  roles.  For  example,  the 
transducers  which  monitor  and  control  vibrations  in  an  aircraft  fuselage  may  also  be  required  to  act  as 
inputs  and  sensors  for  health  monitoring  or  nondestructive  evaluation  of  the  structure.  The  limitations 
on  the  mass  and  size  of  transducers  are  often  relaxed  in  industrial  applications  but  output  requirements 
may  be  more  stringent.  For  example,  the  actuators  employed  for  controlling  vibrations  in  the  cutting 
head  of  a  milling  machine  are  required  to  generate  massive  forces  throughout  a  wide  frequency  range. 

Actuators  and  sensors  comprised  of  smart  or  active  materials  can  meet  these  criteria.  Like  the 
definition  for  smart  structures,  the  definitions  of  smart  or  active  materials  can  vary  somewhat  between 
fields.  We  define  active  actuator  materials  as  those  which  can  convert  electrical,  magnetic  or  thermal 
energy  to  mechanical  energy  while  sensor  effects  are  provided  by  the  opposite  conversion  of  energy.  In 
this  discussion,  we  focus  on  piezoelectric,  electrostrictive,  magnetostrictive  materials  and  shape  memory 
alloys.  The  first  three  compounds  can  be  used  for  both  sensing  and  actuating  while  shape  memory  alloys 
are  used  only  in  the  latter  capacity. 

Piezoelectric  compounds  have  been  the  most  widely  used  transducer  for  smart  structure  design  due 
to  the  fact  that  they  are  lightweight  and  compact,  relatively  inexpensive,  and  exhibit  fairly  linear  field- 
strain  relations  at  low  drive  levels.  Due  to  the  inherent  direct  and  converse  piezoelectric  effects,  they  can 
also  be  employed  as  both  sensors  and  actuators.  At  higher  drive  levels,  however,  they  exhibit  hysteresis 
which  must  be  accommodated  in  applications  that  require  micropositioning  or  precise  control  design. 

In  certain  applications,  electrostrictive  elements  constructed  from  relaxor  ferroelectric  materials  are 
advantageous  over  piezoceramic  patches  due  to  the  fact  that  they  exhibit  minimal  hysteresis  when 
employed  in  the  diffuse  transition  region  near  the  material’s  bulk  Curie  temperature.  This  makes 
them  advantageous  in  systems  including  sonar  transducers  and  deformable  mirrors.  Unlike  piezoelectric 
materials,  electrostrictive  compounds  are  not  poled  and  hence  exhibit  few  aging  effects.  As  detailed 
in  Section  2,  however,  the  constitutive  behavior  of  the  materials  is  highly  temperature-dependent  and 
nonlinear  near  saturation.  Both  aspects  must  be  accommodated  when  designing  control  systems  which 
utilize  the  compounds. 

The  magnetic  analogue  of  the  electrostrictive  compounds  are  magnetostrictive  materials  which  con¬ 
vert  magnetic  energy  to  mechanical  energy  and  vice-versa.  Due  to  the  circuits  required  to  generate 
the  driving  magnetic  fields,  transducers  which  utilize  magnetostrictive  cores  are  currently  larger  and 
more  massive  than  piezoelectric  or  electrostrictive  patches.  The  giant  forces  and  strains  generated 
by  the  transducers,  however,  make  them  advantageous  in  certain  industrial  systems,  and  as  material 
properties  and  transducer  designs  are  refined,  the  scope  of  their  application  should  rapidly  increase. 
From  a  modeling  and  control  perspective,  the  nonlinearities  and  hysteresis  inherent  to  the  materials  at 
moderate  to  high  drive  levels  must  be  accommodated  before  the  materials  can  be  utilized  to  their  full 
potential. 
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Shape  memory  alloys  produce  strains  on  the  order  of  10%  in  response  to  thermal  inputs.  This  makes 
them  candidates  as  actuators  in  applications,  such  as  changing  the  camber  of  a  helicopter  blade,  which 
require  large  strains  or  displacements.  The  modeling  of  these  materials  is  difficult,  however,  due  to 
both  the  large  strains  involved  and  the  inherently  nonlinear  constitutive  behavior  of  the  materials.  As 
detailed  in  Section  4,  this  has  promoted  the  development  of  both  energy-based  and  phenomenological 
models  which  can  be  inverted  for  linear  control  design. 

The  design  of  smart  structures  which  utilize  these  materials  requires  both  the  characterization  of 
their  constitutive  properties  and  the  development  of  coupled  models  which  quantify  their  interaction 
with  underlying  systems.  Control  laws  must  be  compatible  with  the  properties  of  the  sensors  and 
actuators  as  well  as  the  mechanisms  through  which  they  interact  with  the  system.  For  example,  the 
previously  mentioned  actuators  yield  unbounded  (discontinuous)  input  operators  in  the  mathematical 
formulation  of  the  control  problem.  The  extension  of  control  theories  to  this  regime  has  been  been 
completed  in  certain  applications  but  is  lacking  in  general.  Moreover,  all  of  the  active  materials  exhibit 
nonlinear  dynamics  and  hysteresis  at  high  drive  levels.  This  must  be  incorporated  in  both  the  models 
and  control  methods  before  the  materials  can  be  utilized  to  their  full  capability  in  smart  structure  design. 
While  aspects  of  the  nonlinear  design  have  been  addressed,  the  state  of  the  theory  lags  far  behind  that 
for  the  linear  case. 

Because  the  dynamics  of  a  smart  structure  are  dependent  upon  the  attributes  of  the  constituent 
active  materials,  it  is  necessary  to  consider  the  development  of  linear,  nonlinear  and  hysteretic  con¬ 
stitutive  relations,  their  incorporation  in  coupled  system  models,  and  subsequent  control  design  in  a 
concerted  manner.  This  facilitates  the  development  of  models  amenable  to  control  design,  and  the 
formulation  of  control  laws  which  are  compatible  with  the  physical  attributes  of  the  active  sensor  and 
actuator  materials.  Subsequent  discussion  of  the  active  materials  and  their  utilization  in  smart  struc¬ 
tures  follows  this  philosophy.  For  each  compound,  aspects  of  the  material  behavior  which  are  pertinent 
to  model  development  are  discussed  and  appropriate  constitutive  relations  are  presented.  Finally,  a 
variety  of  control  techniques  compatible  with  the  linear  and  nonlinear  models  are  outlined.  While  this 
provides  a  synopsis  of  only  certain  facets  of  model  development  and  control  design  for  smart  structures, 
it  illustrates  issues  which  must  be  addressed  in  order  to  utilize  the  full  potential  of  these  systems. 

2  Piezoelectric  and  Electrostrictive  Materials 

Smart  structure  technology  utilizing  piezoelectric  and  electrostrictive  materials  has  evolved  from  solely 
passive  devices  for  reducing  structural  vibrations  to  highly  integrated  systems  which  provide  the  ca¬ 
pability  for  both  sensing  and  actuating,  adaptively  modifying  system  dynamics  and  monitoring  the 
health  of  the  system.  To  illustrate  facets  of  this  technology  and  the  control  issues  which  it  has  raised, 
we  first  outline  some  prototypical  applications  employing  piezoelectric  and  electrostrictive  transducers. 
A  crucial  aspect  of  control  design  is  the  development  of  appropriate  models  and  this  forms  the  focus 
of  much  of  this  section.  Representative  examples  of  current  control  designs  and  directions  for  future 
research  are  then  summarized  in  Section  2.5. 

2.1  Piezoelectric  and  Electrostrictive  Applications 

2.1.1  Structural  Applications 

Many  of  the  initial  applications  utilizing  both  the  sensing  and  actuating  capabilities  of  piezoelectric  ma¬ 
terials  focussed  on  vibration  attenuation  or  the  reduction  of  structure-borne  noise  in  structural  acoustic 
systems.  To  illustrate  the  first  application,  consider  the  thin  beam  with  surface-mounted  piezoceramic 
patches  depicted  in  Figure  1.  Through  the  application  of  diametrically  out-of-phase  voltages  to  the 
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Figure  1:  (a)  Thin  beam  with  surface- mounted  piezoceramic  patches,  (b)  Bending  moments  generated 
by  out-of-phase  voltages  to  the  patches. 

patches,  bending  moments  are  created  in  the  beam  which  produce  transverse  displacements.  This  pro¬ 
vides  the  actuator  capabilities  for  the  configuration.  The  sensor  capabilities  are  provided  by  the  direct 
piezoelectric  effect  in  which  stresses  in  the  patch  produce  charges  and  corresponding  voltages.  A  truly 
integrated  smart  structure  utilizes  either  self-sensing  actuators  [27]  or  sensor-actuator  pairs  to  provide 
the  observations  and  inputs  required  for  vibration  attenuation. 

While  all  piezoelectric  materials  exhibit  a  certain  degree  of  hysteresis,  it  is  minimal  at  low  drive  levels 
and  can  in  certain  applications  be  reduced  by  charge  or  current  control  [58].  Moreover,  the  feedback 
mechanisms  employed  in  many  control  applications  can  further  reduce  the  hysteresis  exhibited  by  the 
materials.  Hence  the  majority  of  structural  control  applications  utilizing  piezoceramic  transducers  have 
relied  on  linear  models  and  linear  control  methods.  The  assumption  of  linearity  in  piezoelectric  materials 
fails,  however,  at  high  drive  levels  or  high  accuracy  applications  such  as  micropositioning  [34,  35].  In 
such  regimes,  hysteresis  must  be  accommodated  in  both  the  models  and  control  designs.  In  Section  2.3 
we  outline  two  hysteresis  models  which  are  amenable  to  inversion  and  hence  facilitate  the  construction 
of  an  inverse  compensator  for  linear  control  design. 

Relaxor  ferroelectric  materials  such  as  the  electrostrictive  PMN-PT-BT  provide  the  capability  for 
generating  large  strains  with  minimal  hysteresis  if  employed  in  the  diffuse  transition  region.  As  detailed 
in  Section  2.2,  however,  the  relation  between  input  fields  and  generated  strains  is  nonlinear  at  both  low 
drive  levels  and  near  saturation.  While  this  must  be  accommodated  in  control  design,  it  also  provides 
the  possibility  for  designing  adaptive  structures  in  which  coupling  constants  and  sensitivities  can  be 
tuned  by  prescribing  a  bias  field  commensurate  with  a  specific  region  of  the  nonlinear  curve  [63]. 

The  dual  sensor  and  actuator  capabilities  of  piezoelectric  and  electrostrictive  materials  also  provide 
the  possibility  for  health  monitoring  in  smart  structures.  As  detailed  in  [9],  this  can  be  accomplished  by 
driving  the  structure  using  the  actuator  facility  of  the  material  and  sensing  its  response  using  the  direct 
piezoelectric  or  electrostrictive  effects.  The  status  of  the  structure  can  then  be  evaluated  using  various 
criteria.  The  simplest  strategy  is  to  compare  the  response  to  baseline  data  for  the  original  structure. 
If  significant  deviations  are  detected,  more  detailed  analysis  can  be  performed.  If  sufficiently  accurate 
models  are  employed,  this  can  include  the  identification  of  defects  through  least  squares  fits  to  the  data. 
Details  regarding  these  procedures  are  provided  in  [9]  and  included  references. 

In  addition  to  vibration  control  and  health  monitoring,  piezoelectric  and  electrostrictive  actuators  are 
increasingly  employed  to  provide  shape  or  structural  modifications  which  enhance  system  performance. 
A  prototypical  example  from  [49]  is  the  flap  assembly  depicted  in  Figure  2.  In  this  design,  the  strains 
produced  by  stacked  piezoelectric  actuators  are  amplified  by  a  rod  and  cusp  assembly  to  provide  the 
stroke  required  to  rotate  flaps.  This  research  is  directed  at  modifying  the  aerodynamic  properties  of 
aircraft  wings  or  helicopter  rotor  blades. 
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Figure  2:  Piezoceramic  stack  used  for  flap  rotation  (after  [49]). 


2.1.2  Structural  Acoustic  Applications 

A  second  class  of  smart  structures  in  which  piezoceramic  actuators  and  sensors  have  played  a  funda¬ 
mental  role  are  those  involved  in  structural  acoustic  systems.  Structure-borne  noise  arises  in  settings 
ranging  from  aircraft  and  automobiles  to  fields  emanating  from  high  voltage  transformers.  In  all  cases, 
the  unwanted  noise  is  generated  by  structural  vibrations  produced  by  an  adjacent  source  (e.g.,  aircraft 
engines,  impinging  flowfields,  vibrating  machinery,  electromagnetic  cores).  Smart  materials  provide  the 
capability  for  reducing  structure-borne  noise  by  modifying  the  structural  dynamics  to  regimes  which 
couple  less  effectively  with  acoustic  fields  as  depicted  in  Figure  3.  This  capability  can  be  enhanced 
through  the  use  of  models  which  predict  sound  power  levels  as  a  function  of  structural  displacements 
and  velocities.  In  this  case,  piezoceramic  patches  can  be  utilized  to  sense  the  structural  strains  and 
produce  bending  moments  in  accordance  with  a  structural  acoustic  model-based  control  law  [19]. 

The  sound  pressure  levels  encounted  in  originally  considered  aircraft,  automotive  and  industrial 
systems  were  typically  linear,  and  linear  PZT  models  were  employed  for  characterization  and  control 
design.  As  noted  in  Section  2.5,  feedforward  methods  are  currently  employed  for  many  structural 
acoustic  applications.  Because  these  techniques  rely  on  superposition  principles,  they  are  inherently 
linear  and  hence  require  either  linear  models  or  linear  filtering  techniques  such  as  inverse  compensation. 

Current  and  future  structural  acoustic  applications  are  beginning  to  focus  on  regimes  which  are 
highly  nonlinear  and  will  require  the  development  of  nonlinear  models  and  control  techniques.  For 
example,  the  sound  pressure  levels  produced  in  a  space  launch  vehicle  payload  fairing  during  liftoff  or 


Figure  3:  Thin  shell  with  surface- mounted  piezoceramic  patches  which  encloses  and  acoustic  cavity. 
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impinging  on  an  aircraft  bay  at  supersonic  speeds  can  exceed  150  dB  [29,  72,  91].  Smart  materials 
employed  in  these  applications  will  quite  likely  be  operating  in  highly  nonlinear  ranges  which  will 
necessitate  the  use  of  nonlinear  constitutive  and  hysteresis  models  of  the  type  illustrated  in  Section  2.3. 

2.1.3  High  Performance  Transducers 

An  area  of  research  which  has  the  potential  for  significantly  affecting  smart  structure  design  is  the 
development  of  high  performance  transducers  which  employ  piezoelectric  or  electrostrictive  materials. 
Because  these  transducers  can  exhibit  nonlinear  dynamics  and  hysteresis,  their  incorporation  will  also 
necessitate  the  use  of  nonlinear  models  of  the  type  outlined  in  Section  2.3  and  control  methods  com¬ 
mensurate  with  these  models. 

The  number  of  new  high  performance  transducers  that  have  recently  been  developed  is  vast  so  we 
restrict  our  attention  to  three  types  and  refer  the  reader  to  [63]  for  general  discussion  concerning  the 
state  of  the  art  in  this  area. 

The  first  are  THUNDER  actuators  which  are  typically  comprised  of  a  piezoceramic  wafer,  a  metallic 
backing  material,  hot  melt  adhesive  layers,  and  optional  metallic  top  layers  [39,  60].  The  manufacturing 
process  produces  a  characteristic  curved  shape  and  provides  the  actuators  with  the  robustness  necessary 
to  withstand  the  high  voltages  required  for  large  force  and  strain  generation.  Because  the  actuators 
operate  at  extremely  high  drive  levels,  the  relation  between  the  input  field  and  generated  strains  exhibits 
significant  hysteresis  and  nonlinear  constitutive  relations  which  must  be  accommodated  in  models  and 
subsequent  control  design. 

A  second  class  of  composite  transducers  are  flextensional  designs  which  utilize  PZT  or  PMN-PT-BT 
drivers  [21,  36].  Prototypical  designs  are  depicted  in  Figure  4,  and  the  reader  is  referred  to  Section  3.1 
and  [24]  for  discussion  about  analogous  transducers  which  employ  magnetostrictive  drivers.  This  class 
of  transducers  was  originally  developed  for  sonar  projection  but  due  to  its  wide  range  of  hydrodynamic 
response,  it  has  recently  been  used  for  applications  including  oil  exploration  and  underwater  imaging. 
For  high  drive  levels  and  electrostrictive  drivers,  nonlinear  models  must  be  employed  to  characterize 
the  drive  dynamics.  A  second  area  of  modeling  and  optimal  design  involves  the  determination  of 
cap  geometries  which  produce  desired  frequency  responses  for  both  sending  and  receiving  hydrostatic 
signals.  Because  the  transducers  have  both  tunable  actuator  and  sensor  attributes,  they  provide  a 
unique  capability  for  certain  adaptive  smart  structure  applications. 

The  final  class  of  actuators  which  we  mention  are  the  inch  worm  actuators  which  utilize  piezoelectric, 
electrostrictive  or  magnetostrictive  materials  as  drive  elements  and  clamps  [24] .  A  hybrid  design  utilizing 
piezoceramic  and  magnetostrictive  components  is  illustrated  in  Figure  21  of  Section  3.1.  Such  designs 
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Figure  4:  Flextensional  piezoceramic  and  electrostrictive  transducers,  (a)  Piezoceramic  cymbal  actuator 
(after  [21])  and  (b)  PMN-  PT-BT  transducer  (after  [36]). 
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provide  the  capability  for  essentially  unlimited  displacements  with  speeds  currently  on  the  order  of 
1  mm/sec.  As  the  designs  and  drive  electronics  improve,  the  utilization  of  such  devices  in  smart 
structure  applications  should  continue  to  increase. 

2.2  Physical  Properties  of  the  Materials 

To  motivate  the  models  currently  employed  for  smart  structures  employing  ferroelectric  and  relaxor 
ferroelectric  materials,  it  is  necessary  to  understand  certain  characteristics  of  the  material  properties. 
For  example,  the  ionic  biases  produced  during  the  poling  of  piezoelectric  materials  are  the  source  of 
bidirectional  strains  and  quasilinear  behavior  at  low  drive  levels  whereas  the  hysteresis  is  a  direct 
consequence  of  the  domain  structure  inherent  to  ferroelectric  materials.  In  this  section,  we  summarize 
the  material  properties  pertinent  to  subsequent  model  development  and  control  design.  We  consider  the 
commonly  employed  piezoelectric  material  Pb(Zr,Ti)03  and  the  electrostrictive  material  Pb(Mg,Nb)03 
which  are  respectively  referred  to  as  PZT  and  PMN.  For  brevity,  we  focus  on  the  ferroelectric  properties 
of  these  materials  and  refer  the  reader  to  [37]  for  details  concerning  the  physical  behavior  of  relaxor 
ferroelectrics  when  employed  in  the  diffuse  transition  region. 

2.2.1  Ferroelectric  Domains 

Currently  employed  PZT  transducer  materials  are  comprised  of  PbTq.^Os  and  PbZrxO:3  compounds 
with  x  chosen  to  optimize  electromechanical  coupling.  To  simplify  this  discussion,  we  focus  solely  on 
PbTi[_a:0:j  and  refer  the  reader  to  [63]  for  a  general  discussion  of  piezoelectric  material  properties. 

As  illustrated  in  Figure  5,  the  material  distorts  from  a  cubic  to  tetragonal  form  as  it  is  cooled  below 
the  Curie  point.  This  is  manifested  in  the  biasing  of  the  Ti4+  ion  toward  an  O2-  pair  to  minimize  the 


Figure  5:  (a)  Polar  form  of  PbTiOs;  (b)  Potential  energy  as  a  function  of  Ti  position  along  the  X3 
axis;  (c)  Change  in  polarization  AP3  =  ^3303  due  to  an  applied  transverse  stress  ay  (d)  Change  in 
polarization  AP3  =  (I3 \  0\  due  to  an  axial  stress  cj\  (after  [63]). 
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Figure  6:  (a)  Surface  charge  due  to  a  spontaneous  polarization,  (b)  twinned  180°  domains  which  form 
to  minimize  electrostatic  energy. 


potential  energy  and  it  produces  a  spontaneous  polarization  within  the  cell.  As  a  direct  consequence  of 
the  formation  of  spontaneous  polarization  is  the  generation  of  a  surface  charge  and  depolarizing  field 
as  depicted  in  Figure  6.  The  minization  of  electrostatic  energy  then  produces  a  twinned  180°  domain 
structure  where  domains  refer  to  regions  where  dipoles  moments  are  highly  aligned  and  domain  walls 
denote  transition  regions  between  domains.  Furthermore,  the  generation  of  stresses  in  the  material  as  it 
cools  produces  90°  domains  which  minimize  strains  and  the  corresponding  elastic  energy.  We  note  that 
while  polarization  models  must  accommodate  the  switching  of  both  90°  and  180°  domains,  the  former 
is  of  primary  interest  when  quantifying  the  strains  generated  by  the  material. 

The  discussion  thus  far  pertains  to  single  crystals  whereas  the  piezoelectric  materials  used  in  smart 
structure  applications  are  polychrystalline  and,  due  to  the  random  orientation  of  grains  and  domains, 
have  no  net  polarization.  To  produce  the  electrical  anisotropy  which  provides  piezoelectric  compounds 
with  sensor  and  actuator  capabilities,  it  is  necessary  to  pole  the  material  by  applying  an  intense  DC  field 
at  elevated  temperatures  slightly  below  the  Curie  point.  This  partially  aligns  domains  which  produces 
a  net  polarization  in  the  material. 

Two  manifestations  of  the  piezoelectric  effect  for  poled  ferroelectric  materials  are  utilized  in  smart 
structure  applications.  The  direct  effect  involves  the  generation  of  a  charge  when  stress  is  applied  to 
the  material  and  provides  the  material  with  its  sensor  capabilities.  The  converse  effect  reflects  the 
strain  generated  by  an  applied  field  and  is  used  for  actuator  applications.  As  illustrated  in  Figure  5,  the 
applications  of  a  stress  a 3  parallel  to  the  dipole  produces  an  enhancement  of  the  spontaneous  polarization 
whereas  the  application  of  a  perpendicular  stress  <ji  decreases  the  spontaneous  polarization.  Hence  the 
piezoelectric  coupling  coefficient  dy.i  is  positive  while  ei-j  1  is  negative. 

2.2.2  Hysteresis 

An  important  manifestation  of  the  ferroelectric  nature  of  piezoelectric  and  electrostrictive  materials  is 
the  inherent  hysteresis  observed  at  moderate  to  high  input  levels.  A  typical  hysteresis  relation  between 
the  input  field  E  and  polarization  P  for  PZT5A  is  plotted  in  Figure  7.  The  nearly  vertical  burst  region 
of  the  curve  is  attributed  to  180°  domain  nucleation  and  growth.  Near  saturation,  the  material  acts  as 
a  single  domain  and  subsequent  changes  in  polarization  are  due  to  the  extension  of  dipoles  through  ion 
displacement. 

Sigmoidal  hysteresis  of  the  form  depicted  in  Figure  7  is  typically  attributed  to  the  restriction  of 
domain  wall  movement  by  inclusions  or  pinning  sites  inherent  to  the  material.  These  sites  can  be 
caused  by  impurities,  second-phase  materials,  or  stress  nonhomogeneities  which  in  many  cases  are  not 
only  unavoidable  but  are  necessary  for  the  material’s  performance.  The  pinning  of  domain  walls  at  these 
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Figure  7:  Hysteresis  observed  in  a  PZT5A  wafer  in  response  to  a  1600  V  input. 

locations  results  from  the  reduction  in  electrostatic  energy  which  occurs  at  the  inclusions.  A  complete 
discussion  of  the  mechanisms  which  underlie  domain  wall  pinning  and  movement  is  beyond  the  scope  of 
this  presentation  and  we  refer  the  reader  to  [80]  for  details  and  references  concerning  these  phenomena. 

For  subsequent  model  development,  we  simply  note  that  the  motion  of  domain  walls  pinned  at 
inclusions  is  comprised  of  two  primary  components:  bending  and  translation.  Bending  occurs  at  low 
input  field  levels  before  domain  walls  have  attained  sufficient  energy  to  move  out  of  local  energy  wells 
[15,  30,  54],  This  motion  is  reversible  and  accounts  for  a  large  percentage  of  the  polarization  change 
in  hard  PZT  materials.  At  higher  field  levels,  the  local  energy  barriers  are  overcome  and  domain  walls 
translate  to  a  remote  pinning  site  [15].  This  provides  an  irreversible  component  to  the  polarization 
which  is  present  in  both  hard  and  soft  materials  but  plays  a  more  major  role  in  soft  materials  where 
lower  energies  are  required  to  rotate  domains.  These  reversible  and  irreversible  components  to  the 
magnetization  will  be  quantified  in  Section  2.3  to  provide  a  hysteresis  model  for  ferroelectric  materials. 

2.3  Material  Models 

To  characterize  the  relations  between  applied  fields  and  the  polarization  which  results  in  ferroelectric  and 
relaxor  ferroelectric  materials,  we  consider  three  classes  of  models.  The  linear  models  are  applicable  for 
piezoelectric  materials  at  low  drive  levels  where  hysteresis  and  saturation  phenomena  are  minimal.  The 
second  class  of  models  includes  higher  order  representations  which  are  suitable  for  relaxor  ferroelectric 
operating  in  the  diffuse  transition  region  where  hysteresis  is  minimal.  Finally,  we  outline  two  methods 
for  quantifying  the  hysteresis  observed  in  piezoelectric  materials  at  high  drive  levels  and  electrostrictive 
materials  operating  in  the  ferroelectric  range.  The  first  hysteresis  model  is  based  on  electrostrictive 
domain  principles  while  the  second  is  a  Preisach  approach  which  is  phenomenological  in  nature.  Both 
techniques  are  amenable  to  inversion  which  facilitates  linear  control  design  through  the  construction  of 
an  inverse  compensator. 

2.3.1  Linear  Model 

At  low  input  field  levels,  piezoelectric  materials  exhibit  bidirectional  strains  and  nearly  linear  relations 
between  the  field,  voltage,  or  polarization  and  the  resulting  strain.  Similarly,  the  electric  fields  generated 
by  the  direct  piezoelectric  effect  exhibit  a  nearly  linear  dependence  upon  the  applied  stress.  The 
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development  of  linear  models  relating  the  dielectric  and  elastic  behavior  of  piezoelectric  materials  is 
attributed  to  Voigt,  and  the  resulting  coupled  models  are  used  for  smart  material  characterization  in 
numerous  applications  which  require  low  input  fields  or  voltages. 

As  detailed  in  [61,  64],  the  assumption  of  linear  dielectric  and  elastic  behavior  respectively  yields 
the  relations 

A  =  dijkOjk  +  eijEj 

Cjj  4“  dkijEk 

as  models  for  the  direct  and  converse  piezoelectric  effects.  Note  that  the  Einstein  summation  convention 
is  employed  in  both  relations  so  that  terms  with  repeated  indices  are  summed.  Here  D  and  E  denote 
the  usual  electric  fields,  a  and  e  respectively  denote  the  applied  stress  and  resulting  strain,  e"-  denotes 
components  of  the  dielectric  permittivity  tensor  measured  at  constant  stress,  and  sE  is  the  elastic 
compliance  tensor  measured  at  constant  electric  field.  The  first  term  on  right  side  of  the  converse 
relation  is  due  to  Hooke’s  law  while  the  second  describes  the  strains  which  occur  in  a  material  which 
is  free  to  deform  in  the  absence  of  stresses.  Finally,  we  note  that  the  relation  Dt  =  e^Ej  in  the  direct 
model  should  be  used  with  caution  for  ferroelectric  materials  since  the  permittivity  is  not  only  nonlinear 
but  is  a  multi-valued  map  at  high  input  levels. 

The  tensor  formulations  in  (1)  provide  a  complete  description  of  the  material  behavior  for  general 
piezoelectric  compounds.  The  number  of  coefficients  can  be  reduced  significantly,  however,  by  invoking 
elastic  symmetries  and  electric  symmetries  due  to  poling.  When  combined  with  a  change  of  index,  this 
provides  a  matrix  system  suitable  for  typical  smart  material  applications. 

We  first  note  that  for  linear  elastic  materials,  the  stress  and  strain  tensors  are  symmetric  so  that 
Oij  =  Uji  and  etJ  =  eji.  The  symmetry  in  stress  and  strain  implies  that  d^k  is  symmetric  in  j  and  k 
thus  leaving  18  independent  coefficients. 

To  formulate  the  system  as  a  matrix  equation,  the  indices  jk  are  replaced  by  a  single  index  m 
according  to  the  following  convention: 

11  22  33  23,  32  31,  13  12,  21  Tensor 

1  2  3  4  5  6  Matrix 

The  coefficients  dijk  are  then  written  as  d.l;jk  =  drrn  when  m  =  1,2,3  and  dijk  —  \dim  when  m  =  4, 5, 6. 
The  factor  of  1  /2  accommodates  the  symmetry  in  j  and  k.  A  similar  convention  is  used  to  formulate 
the  stress,  strain  and  compliance  tensors  as  matrices  (e.g.,  [23,  71]). 

The  number  of  coefficients  dim  and  stress-strain  elements  of  interest  are  further  reduced  when  we 
consider  a  material  which  is  poled  in  the  x%  direction.  As  noted  in  [23],  conical  symmetry  dictates  that 
in  this  case,  all  of  the  piezoelectric  coefficients  are  zero  except  d%\  =  Vi 2 ,  Vi 3  and  di§  =  dn  . 

The  utilization  of  the  structural  and  electrical  symmetries  along  with  matrix  re-indexing  then  yields 
the  system 
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summarizing  the  converse  and  direct  piezoelectric  effects.  In  matrix  form,  the  components  can  be 
written  as 


D  =  da  +  eaE 

(2) 

e  =  sEa  +  <F E  . 

Finally,  by  utilizing  the  relation  Pj  =  toXijEj ,  where  Xij  is  the  dielectric  susceptibility  tensor  and 
eo  is  the  permittivity  of  a  vacuum,  one  can  also  express  the  strain  in  terms  of  the  polarization  through 
the  relation 


eij  —  sfjk£ake  +  dejkttfPk  ■  (3) 

This  formulation  proves  advantageous  when  modeling  ferroelectric  hysteresis  in  the  relation  between 
the  field  and  polarization. 

2.3.2  Higher  Order  Anhysteretic  Models 

The  linear  models  (l)-(3)  are  valid  only  for  poled  materials  with  small  to  moderate  applied  stresses 
or  fields.  For  general  relaxor  ferroelectrics  and  higher  input  levels,  quadratic  electrostrictive  effects, 
saturation  phenomena  and  hysteresis  are  significant  and  must  be  incorporated  in  the  models.  We 
summarize  here  higher  order  models  for  hysteresis-free  (anhysteretic)  behavior  and  address  hysteresis 
in  the  next  section. 

When  formulated  in  terms  of  the  polarization,  the  converse  and  direct  effects  in  centrosymmetric 
materials  (e.g.,  electrostrictive  materials)  can  be  modeled  by  the  relations 


Ei  —  2QkiijO'kiPj  +  P(Pj,T) 

E  (4) 

&ij  —  ^ijki^kl  T  Qijk(PkPi  ■ 

Here  Q  is  a  fourth  rank  tensor  of  electrostrictive  coefficients  while  P(p,  T)  is  a  nonlinear  function  of 
the  polarization  and  temperature  T  which  incorporates  the  low  field  quadratic  behavior  as  well  as  the 
saturation  effects  which  occur  at  high  drive  levels. 

To  illustrate  the  derivation  of  appropriate  functions  P(P,  T)  and  to  provide  the  framework  re¬ 
quired  for  subsequent  hysteresis  models,  we  consider  the  modeling  of  hysteresis  free  (anhysteretic) 
PMN  through  the  application  of  Boltzmann  statistics.  As  detailed  in  [79,  80],  the  potential  energy  for 
a  dipole  moment  p0  in  an  electric  field  E  is 

£  =  — po  •  E  =  —poE  cos  9 

where  po  —  |po|,  E  —  |E|.  The  balancing  of  thermal  and  electrostatic  energy  then  yields  the  probability 

H(£)  =  Ce~£lkBT 


that  a  dipole  occupies  the  energy  state  £.  Here  k u  denotes  Boltzmann’s  constant  and  C  is  a  constant 
which  is  specified  in  a  manner  which  ensures  that  integration  over  all  possible  configurations  yields  the 
total  number  of  moments  per  unit  volume.  The  assumption  that  the  orientation  of  cells  can  be  in  any 
direction  yields  the  classical  Langevin  relation 


P  —  P 

J  an  —  1  s 


coth  (  — 
a 


(5) 
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where  Ps  is  the  saturation  polarization  and 


Ee  =  E  +  aP  +  2Q33P  g  (6) 

denotes  the  effective  field  acting  at  the  domain  level.  For  constant  temperature,  the  parameters  a  and 
a  are  given  by  a  =  E/Ps  and  a  =  ET  /  (3TC)  where  E  denotes  a  scaling  electric  field  and  Tc  is  the  Curie 
temperature.  A  second  model  which  results  from  the  assumption  that  dipoles  align  only  in  the  direction 
of  the  field  or  opposite  to  it  is  the  Ising  spin  relation 

Pan  =  -Pstanh  (V 

where  a  =  ET/TC.  The  behavior  of  the  two  models  is  illustrated  in  Figure  8.  We  also  refer  the  reader 
to  [67]  for  an  alternative  approach  to  modeling  the  anhysteretic  polarization. 

To  obtain  a  constitutive  relation  commensurate  with  (4),  the  characterization  must  be  polarization 
based.  The  necessary  inversion  can  be  accomplished  for  the  Ising  spin  model  which  yields 


ET  P 

E  =  —2Q^Pa  —  aP  +  — -  arctanh— 

C  *  S 

Hence  for  this  model,  the  function  T  in  (4)  is  given  by 


(8) 


ET  P 

P{T,  P )  =  —aP  H - arctanh—  . 


Further  details  concerning  the  derivation  and  attributes  of  this  model  are  provided  in  [79,  80]. 


2.3.3  Hysteresis  Models  -  Domain  Wall  Model 

As  detailed  in  Section  2.2,  hysteresis  in  the  relation  between  the  applied  field  and  resulting  polarization 
is  an  inherent  property  of  all  ferroelectric  materials.  We  summarize  here  two  techniques  for  modeling 
hysteresis  in  high  drive  level  applications.  The  first  is  based  on  domain  energy  relations  while  the  second 
illustrates  a  phenomenological  Preisach  characterization. 
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Figure  8:  (a)  Ising  spin  and  Langevin  models  for  the  anhysteretic  polarization;  (b)  Fit  to  600  V,  800  V, 
1000  V  and  1600  V  PZT5A  data  using  the  model  (14)  with  one  set  of  parameters. 
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The  domain  model  is  based  on  the  observation  that  due  to  inclusions  in  the  material,  changes 
in  polarization  resulting  from  domain  wall  movement  have  irreversible  and  reversible  components.  To 
quantify  the  irreversible  polarization,  it  is  noted  in  [79,  80]  that  the  polarization  level  for  a  given  effective 
field  can  be  expressed  as  that  for  the  ideal  case  minus  losses  required  to  break  pinning  sites.  This  yields 
the  relation 


Pirr  —  Pan,  k 


dPir \ 
dEe 


(10) 


where  Pan  is  defined  by  (5)  or  (7)  and  k  quantifies  the  average  energy  required  to  break  pinning  sites. 
The  formulation  of  (10)  in  terms  of  the  applied  field  E  yields  the  differential  equation 


dp. 


dE 


(11) 


where  the  parameter  6  —  sign  (dE)  ensures  that  the  energy  required  to  break  pinning  sites  always 
opposes  changes  in  polarization.  The  inclusion  of  the  parameter  6  given  by 


S  = 


1  ,  {dE  >  0  and  P  <  Pan}  or  {dE  <  0  and  P  >  Pan} 
0  ,  otherwise 


(12) 


enforces  that  solely  reversible  polarization  changes  take  place  when  the  field  reversed  at  saturation  until 
the  anhysteretic  polarization  is  reached. 

The  second  component  of  the  polarization  is  the  reversible  polarization  which  models  the  effects  of 
domain  wall  bending.  To  first  approximation,  this  is  modeled  by  the  relation 


Prev  —  c( Pan  Pirr ) 

where  c  is  a  parameter  which  must  be  estimated  for  the  specific  application  (see  [79,  80]). 
The  total  polarization  is  then  given  by 


(13) 


P  —  Prev  +  Pit 


(14) 


To  implement  the  model,  the  effective  field  for  a  given  field  and  irreversible  polarization  level  is  computed 
using  (6).  This  value  is  then  employed  in  either  (5)  or  (7)  to  compute  the  corresponding  anhysteretic 
polarization.  The  subsequent  irreversible  polarization  is  determined  by  numerically  integrating  (11). 
The  total  polarization  is  then  specified  by  (14). 

The  experimental  validation  of  the  model  for  PMN  transducers  operating  in  the  ferroelectric  range 
is  illustrated  in  [79,  80]  while  validation  for  various  PZT  compounds  is  addressed  in  [81,  82].  The 
predictive  capability  of  the  model  is  illustrated  in  Figure  8b.  The  parameters  a,a,  c,  k  and  Ps  were 
estimated  through  a  least  squares  fit  to  600  V  data  and  the  model,  with  the  same  parameters,  was  used 
to  predict  the  polarization  resulting  from  800  V,  1000  V  and  1600  V  inputs.  The  formulation  of  the 
model  in  terms  of  electrostatic  energy  principles  provides  it  with  the  capability  for  accurately  predicting 
the  polarization  throughout  the  operational  range  given  only  specified  input  levels. 

A  second  advantage  of  this  model  lies  in  the  fact  that  it  can  be  inverted  through  consideration  of  a 
complementary  ODE.  This  permits  the  construction  of  an  inverse  compensator  which  can  be  employed 
as  a  filter  before  the  actuator  to  facilitate  linear  control  design.  Details  regarding  the  construction  of 
the  inverse  are  provided  in  Section  3.3  where  the  analogous  magnetic  model  is  outline. 
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Figure  9:  (a)  Preisach  plane  S  and  (b)  Restricted  Preisach  plane  S a- 


2.3.4  Hysteresis  Models:  Preisach  Model 

A  second  technique  for  modeling  hysteresis  in  a  variety  of  smart  material  systems  is  based  on  the  ap¬ 
proximation  of  the  multi-valued  hysteresis  map  by  a  parallel  collection  of  Preisach  kernels.  This  provides 
a  phenomenological  or  empirical  characterization  of  the  input /output  relations  which  circumvents  un¬ 
modeled  or  unknown  physics.  An  advantage  of  this  technique  lies  in  its  generality  as  evidenced  by  the 
fact  that  the  same  model  can  be  used  to  characterize  the  hysteresis  in  piezoelectric  materials  or  shape 
memory  alloys  (see  Section  4.1).  A  second  advantage  is  the  fact  that  the  resulting  Preisach  operator 
can  be  inverted  to  facilitate  linear  control  design.  A  disadvantage  of  the  approach  is  that  it  typically 
requires  the  estimation  of  a  large  number  of  nonphysical  parameters  and  it  does  not  utilize  or  provide 
physical  insights  regarding  the  material  dynamics  of  the  system. 

We  focus  on  the  Krasnolselskii-Pokrovskii  (KP)  operator,  also  termed  a  generalized  Preisach  oper¬ 
ator,  since  it  provides  continuity  with  respect  to  both  time  and  parameters  and  is  amenable  to  approx¬ 
imation  [6,  7].  The  reader  is  referred  to  [34,  35]  for  a  discussion  of  adaptive  control  design  for  classical 
Preisach  operators  defined  in  terms  of  piecewise  constant  kernels. 

General  Preisach  operators  are  defined  by 

[pu(v,ZW)  =  [_[ks{v,Z{s))]{t)dn{s)  (15) 

Js 


where  fi  €  M  is  a  finite  signed  Borel  measure  on  the  closed  Preisach  plane  S  =  {(si,S2)  |  <  $2}  as 

depicted  in  Figure  9.  The  input  is  denoted  by  v  while  denotes  initial  states  of  the  system.  In  classical 
Preisach  formulations,  the  kernel  ks  is  defined  through  a  linear  combination  of  piecewise  constants 
whereas  continuous,  piecewise  linear  functions  are  employed  in  the  KP  kernel  to  ensure  continuity  with 
respect  to  time  and  parameters.  The  characterization  of  hysteresis  using  either  kernel  requires  the 
identification  of  an  appropriate  measure  /i(.s)  which  acts  as  a  weight  on  the  kernel. 

To  define  the  KP  kernel,  we  consider  translates  rSl  =  r(v  —  s  j )  and  rS2  =  r(v  —  S2)  of  a  Lipschitz 
continuous  ridge  function  r(v)  as  depicted  in  Figure  10.  For  time  intervals  [1/..  1 ,  t^]  where  the  input  v 
is  monotone,  a  monotone  operator  is  recursively  defined  by 


[77(u,77A;_i)](t) 


max{77.£-i,  r(v(t)  —  S2)}  if  v  is  non-decreasing 
min{77yfc_i, r(v(t)  —  si)}  if  v  is  non-increasing 


where 


Kk 


K{v,Kk-i){tk)  ,  k  =  2 

770  =  C,  k  =  1,|  G  {-1,1} 
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(a) 


Figure  10:  (a)  Ridge  function  r(v  —  s).  (b)  Hysteresis  envelop  provided  by  the  translates  r\  and  r-i  of 
the  ridge  function  r(v). 


defines  the  values  of  7 Z  at  times  f*.  The  KP  kernel  is  then  defined  recursively  on  each  subinterval  by 

[ks{v,$,)](t)  =  [niv.TZk-i)]^)  ,  f  6  [4-1,4]-  (16) 

A  typical  path  for  ks  is  depicted  in  Figure  11  while  details  concerning  the  properties  of  the  kernel  are 
provided  in  [6,  7]. 

To  numerically  implement  the  model,  the  measure  n(s)  is  approximated  by  the  expansion 

m 

l^m  =  ani^ni  (17) 

»:  1 

where,  as  depicted  in  Figure  9,  n*  are  nodes  in  the  restricted  Preisach  plane  S a  =  {(si,  sa)  | smin  < 
si  <  S2  <  smax},  «ni  are  weights,  and  6rii  denotes  the  Dirac  measure  having  an  atom  at  n,.  The  KP 
operator  P/(  is  then  approximated  by  Pm  given  by 

m 

[-Pm(^,0](i)  =  ■  (18) 

i—  1 

As  detailed  in  [6,  7,  32],  the  continuity  and  convergence  properties  of  the  KP  operator  lead  to  well- 
posedness  of  the  corresponding  identification  problem  and  provide  a  framework  which  is  amenable  to 
implementation. 

An  attribute  of  both  the  classical  Preisach  and  KP  operators,  which  is  crucial  to  control  design, 
is  the  fact  that  they  can  be  inverted.  We  summarize  here  details  regarding  the  inversion  of  the  KP 
operator  and  refer  the  reader  to  Section  1.4  for  an  example  illustrating  its  use  in  control  design. 


Figure  11:  (a)  Piecewise  monotone  input  v(t).  (b)  Output  from  the  KP  kernel  in  response  to  the  input. 
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Figure  12:  (a)  Plant  with  input  actuator  hysteresis,  (b)  Approximate  inverse  compensator  P,m  1  for  the 
hysteresis. 

The  subsequent  development  of  the  inverse  operator  P,m  1  follows  that  in  [32],  We  first  note  that  if 
Ud(t)  denotes  the  desired  actuator  input,  the  process  of  inverse  compensation  can  be  represented  by 

u  =  =  Pm{Pm  =  ud 

as  illustrated  in  Figure  12. 

With  the  definition  A  Ud  =  u(i{t)  —  u,i(to).  it  then  follows  that  changes  in  the  desired  actuator  input 
on  the  interval  [to,  ^i]  are  given  by 

Aitd  =  Pm{y ,£)  Pm[y  o,£) 

m 

=  ^  ~  [kniiv0i  Olni 

i=  1 

=  Y  "  n«)  "  r(w  “  n*)]  am  ■ 

nies  a 

For  the  piecewise  linear  ridge  function 

{1  .  .s2  <  v  a 

—  1  +  2(v  —  s2)/a  ,  v  —  a  <  s2  <  v 

-1  ,  v  <  s2 

depicted  in  Figure  10,  the  change  in  the  kernel  at  each  node  s2  =  nt  is  given  by 

0  ,  rii  <  vo  —  a 

r{v  -  rii)  -  r(v0  -m)  =  l  A rUi  ,  v0  -  a  <  n*  <  v  (19) 

0  ,  v  <rii. 

It  then  follows  that 

^ ud =  ani  —  ani  +  yy  Arn!.ani 

n,e5+  nj£5A 

where  _ 

S~  =  {rii  E  S a  | r(v  -  rii)  =  -1} 

S+  =  { rii  E  Sa  | r(v  -  rii)  =  1} 

SA  =  { rii  £  Sa  |  —  1  <  r(v  —  rii )  <  1}  • 

If  A u  is  less  than  or  equal  to  the  distance  between  adjacent  nodes  n*,  then 

A 7-j  =  2(v  -  vq)  ja 
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and  (19)  can  be  solved  for  v  to  obtain  the  relation 

Pm1^^!)  =  v0  +  a(Aud-  ^2  am  +  Y  a^)/(2  Y  a^) 

n,es+  UiES -  n,e5A 

for  the  inverse  operator.  As  detailed  in  [32],  an  iterative  procedure  provides  the  inverse  if  A u  is  greater 
than  the  distance  between  nodes.  Experimental  validation  of  this  modeling  and  inversion  procedure  for 
quantifying  hysteresis  for  a  piezoceramic  stack  is  reported  in  [32], 

2.4  System  Models 

The  discussion  in  Section  2.1  indicates  certain  smart  structure  applications  which  utilize  piezoceramic 
and  electrostrictive  actuators  while  the  models  of  Section  2.3  characterize  certain  material  properties 
for  these  actuator-sensor  compounds.  In  this  section,  we  illustrate  the  incorporation  of  the  constitutive 
relations  from  Section  2.3  in  prototypical  structural  and  structural  acoustic  systems,  and  the  formulation 
of  the  system  model  in  a  format  suitable  for  control  implementation.  For  brevity,  we  consider  a  thin 
beam  model  and  a  coupled  2-D  model  as  prototypes  for  structural  and  structural  acoustic  systems. 
While  both  represent  simplified  scenarios,  they  provide  the  framework  required  for  modeling  more 
complex  systems  as  a  prelude  to  control  design. 

2.4.1  Beam  with  Piezoceramic  Actuators  and  Sensors 

Consider  a  cantilevered  thin  beam  with  a  pair  of  surface-mounted  piezoceramic  actuators  as  depicted  in 
Figure  13.  For  modeling  purposes  the  beam  is  assumed  to  have  length  l.  width  6,  and  thickness  h.  The 
density,  Young’s  modulus  and  Kelvin- Voigt  damping  coefficient  for  the  beam  and  piezoceramic  patches 
are  respectively  denoted  by  p&,.S&,cd6  and  ppe.  E,pc.  c/j)je .  The  transverse  beam  displacement  is  given 
by  w  while  p(t.  x )  denotes  an  exogenous  surface  force  to  the  beam.  We  note  that  for  this  configuration, 
the  linear  piezoelectric  equations  (2)  reduce  to 

D:i  -  d31cr3  +  e3E3  Direct 

1  ,  ,  V  (20) 

e3  =  —  cr3  +  «3i  j—  Converse 

J2/pe  ripe 

where  V  denotes  the  voltage  to  the  patches  and  we  have  employed  the  relation  E  —  V/hpe  in  the 
converse  relation. 


Figure  13:  Thin  beam  with  surface-mounted  piezoceramic  patches. 
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As  detailed  in  [9] ,  force  and  moment  balancing  yields  the  strong  form  of  the  Euler-Bernoulli  equations 

Pfa)^fa-fa  *)  +  9  fa?*  &  x">  =  Pfa  +  9  fa  x) 

wfa°)  =  ^fa°)  =  Q  (21) 

=  ^p(t,€)  =  0 

along  with  appropriate  initial  conditions,  as  a  model  for  the  transverse  beam  dynamics.  The  composite 
density  and  internal  bending  moment  are  given  by 

Pfa)  =  Pbhb  +  2  PpebhpeXpztfa) 

d2w  d3w  ^  ^ 

Mint  fax)  =  El(x)-^fax)  +  cDI-^^fax) 

where  the  characteristic  function  Xpzt  delineates  the  location  of  the  rods  and 

,  T ,  ,  Efa^b  26 
El  fa)  —  —  V  —EpedsXpztfa) 

12  ,  d  (23) 

T,  .  CDhh3b  26 
CDlfa)  =  12  +  yCDpe03  Xpfax:) 

with  «3  =  (6/2  +  hpe )3  —  63/8.  The  first  term  in  the  definition  of  El  is  the  usual  stiffness  parameter 
for  a  uniform  beam  while  the  second  is  due  to  the  passive  material  contributions 

rh/2+hPe  26 

2  /  Epee3bzdz  =  —Epea3 

Jh/2  -J 

from  the  converse  patch  relation.  The  factor  of  two  results  from  invoking  the  symmetry  of  the  two 
patches.  Integration  of  the  strain  rate  for  the  beam  and  patch  material  yields  the  Kelvin- Voigt  damping 
parameter. 

For  the  case  when  the  patches  are  driven  diametrically  out-of-phase,  the  integration  of  the  free  strain 
component  of  the  converse  relation  yields  the  external  moment 

r  h/2-\-hpe 


rri/  zfaripe 

=  2  /  d>3iEpeb(V / hpe) z  dz 

Jh/2 


=  KB  XpztV  fa) 

where  KB  =  Epebfai(h  +  hpe)  and  V(t)  is  the  voltage  to  the  patches.  In  applications,  the  actuator 
coefficient  Kb  is  typically  identified  for  a  given  structure  through  a  least  squares  fit  to  data. 

To  obtain  a  weak  form  of  the  model,  we  take  the  state  to  be  the  displacement  w  in  the  state  space 
X  =  L2(0,  l)  with  the  inner  product 

fafa)x  =  [  pHdx. 

Jo 

The  space  of  test  functions  is  taken  to  be  V  =  H2(0,£)  =  {<j>  E  H2(0,£)  \  (f)(0)  =  cf)'( 0)  =  0}  with  the 
inner  product 

fa 

fafa)v  =  /  El  fa  fa"  dx  . 
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It  should  be  noted  that  with  these  choices,  V  is  continuously  and  densely  embedded  in  X.  Hence  one 
has  the  Gelfand  triple  V  <->■  X  ~  X*  <->  V*  with  the  pivot  space  X. 

A  weak  form  of  the  model  is  then  given  by 


nt  r£  rl  r£ 

/  pwtpdx+  /  Mint4>"  dx  =  /  Mexti>"  dx  +  /  pipdx  (25) 

Jo  Jo  Jo  Jo 

for  all  if)  e  V. 

The  equations  (21)  and  (25)  with  moments  given  by  (22)  and  (24)  model  the  actuator  contributions 
to  the  structure.  The  sensor  attributes  are  modeled  using  the  direct  piezoelectric  relation  in  (20).  For 
a  beam  in  which  strains  are  solely  due  to  bending,  the  stress  at  a  point  (x,zc),  where  zc  =  (h  +  hpe)/2 
is  the  distance  from  the  neutral  axis  of  the  beam  to  the  center  of  the  patch,  is  approximated  by 

.  d2w 

CT3  (t,x)  =zcEpe—I. 

In  the  absence  of  an  applied  field  ( E3  =  0),  the  direct  piezoelectric  relation  then  yields 

d2w 

D3  =  d3izcEpe-^--j  .  (26) 

The  charge  q(t)  developed  on  the  patch  as  the  material  strains  is  computed  by  integrating  (26)  over  the 
area  of  the  patch  to  obtain 


fX2  d2w 

q(t)  =  b  d3izcEpe—. Tdx 

J  X\ 


=  Epgdsibz, 


pe  dx2 
dw 


dw 


Finally,  the  voltage  Vp  generated  by  the  patch  is  computed  by  dividing  the  charge  by  the  capacitance 
Cpe  to  obtain 


Vp(t)  =  Ks 


dw  dw 


(27) 


where 

j£S  _  Eped3ib(h  +  hpe) 

2  Cpe 

Hence  the  patch  can  be  used  to  sense  accumulated  strain. 

The  previous  discussion  illustrates  the  use  of  the  direct  and  converse  piezoelectric  effects  for  em¬ 
ploying  piezoceramic  patches  as  sensors  and  actuators  in  smart  structures.  Additionally,  the  dual 
piezoelectric  properties  provide  the  patches  with  the  capability  for  simultaneous  sensing  and  actuating. 
Self-sensing  piezoceramic  actuators  were  developed  somewhat  concurrently  by  Dosch  et  al.  [27]  and 
Anderson  and  Hagood  [2]  while  adaptive  sensoriactuators  were  developed  by  Cole  and  Clark  to  address 
difficulties  associated  with  drifting  capacitance  [20] .  Whether  the  patches  are  employed  as  isolated  sen¬ 
sors  and  actuators  or  as  collocated  sensoriactuators,  they  provide  significant  potential  as  transducers  in 
smart  structures. 


18 


2.4.2  Finite  Dimensional  Model 


It  is  necessary  to  approximate  the  dynamics  of  the  infinite-dimensional  model  (25)  before  it  can  be 
employed  for  either  simulation  purposes  or  control  design.  We  employ  a  Galerkin  approximation  in  the 
spatial  variable  to  obtain  a  semidiscrete  ODE  system  in  time  which  is  amenable  to  control  formulation. 
Specifically,  the  spatial  basis  is  taken  to  be  where  4>j{x)  denotes  the  jth  cubic  U-spline  modified 

to  satisfy  the  fixed  left  boundary  condition  [68].  Approximate  solutions 

m+ 1 

Wm{t,x)  =  WjiWjix) 

3= 1 

are  then  considered  in  the  subpace  Vm  =  span{<^}.  To  obtain  a  vector  ODE  system,  the  infinite 
dimensional  system  (25)  is  restricted  to  Vm  and  posed  in  first-order  form  to  yield 

x(t)  =  Ax(t)  +  Bu(t )  +  g(t) 
z(0)  =  Xq  . 


(29) 


(28) 


The  state,  control  input,  and  control  matrices  are  respectively  defined  by  x(t)  =  [u>i(f),  ■  ■  ■  ,  wm+i(t). 
wi(t),  ■  ■  ■  ,wm+i(t)],  u(t)  =  [tq(f),  ■  ■  ■  ,  us(f)]  and  B  =  KDb.  The  component  system  matrices  have  the 
form 


A  = 


0  I 

Q~XK  Q~xCd 


b  = 


0 

Q~lb 


(30) 


where 


9{t) 


0 

Q_1p(t) 


[Q\ij  ~ 

f 

/  p4>i4>j  dx 

"e 

% 

II 

(S 

<0 

J  mag 

ft 

[K\rj  = 

/  EBp'Icj)1!:  dx 

[p{t)]i=  /  p{t,x)4>idx 

Jo 

Jo 

r(- 

[' CD\ij  = 

/  cDI^((/;  dx . 

Jo 

Details  regarding  the  construction  of  the  finite  dimensional  system  can  be  found  in  [9,  77]. 


2.4.3  Nonlinear  Actuator  Models 

The  linear  system  model  (29)  is  based  on  the  assumption  that  the  linear  constitutive  relations  (20) 
adequately  model  the  actuator  and  sensor  effects  of  the  patches.  As  detailed  in  Section  2.2,  however,  both 
piezoceramic  and  electrostrictive  materials  exhibit  nonlinear  constitutive  behavior  and  varying  degrees 
of  hysteresis  at  high  drive  levels  which  necessitates  consideration  of  the  hysteresis  models  discussed  in 
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Section  2.3.  The  incorporation  of  domain  wall  or  Preisach  hysteresis  models  with  the  beam  equations 
yields  the  system 

x{t)  =  Ax(t)  +  [5(it)](f)  +g{t) 

x(0)  =  Xq 

upon  approximation.  The  nonlinear  input  [i?(ri)](t)  =  [P(u)](t)b  where  6,  given  by  (30),  incorporates 
the  spatial  location  of  the  actuator  and  [P(u)](t)  quantifies  the  hysteresis  and  nonlinear  behavior. 
Details  regarding  the  construction  of  B(u )  for  the  analogous  magnetic  domain  wall  hysteresis  model  are 
provided  in  Section  3.5  while  it  construction  for  the  Preisach  operator  is  illustrated  in  the  context  of 
shape  memory  alloys  in  Section  4.1. 


2.4.4  Structural  Acoustic  System 


The  second  example  which  we  consider  is  a  prototypical  structural  acoustic  system  in  which  piezoceramic 
actuators  mounted  on  the  structure  are  employed  as  actuators  and  sensors.  To  simplify  the  discussion, 
we  consider  a  rectangular  cavity  with  a  vibrating  beam  at  one  end  as  depicted  in  Figure  14.  The 
remaining  cavity  walls  F  are  assumed  to  reflect  all  acoustic  vibrations.  A  force  p  modeling  an  exogenous 
pressure  field  drives  the  system  and  produces  unwanted  noise  in  the  cavity  Fl. 

As  detailed  in  [9],  for  sound  pressure  levels  on  the  order  of  those  typically  encountered  in  many 
aeronautical  and  automotive  applications,  the  acoustic  dynamics  in  a  medium  with  density  pj  can  be 
modeled  by  the  wave  equation 


^■  =  c2A(f)  ,  (i,j)6fl 
V<p  ■  h  —  0  ,  (x,y)  €  T 


(32) 


where  cf>  denotes  the  acoustic  potential,  c  is  the  speed  of  sound,  p  =  Pf-$  is  the  pressure  and  V<^>  =  —u 
is  the  field  velocity.  The  balancing  of  velocity  at  the  beam  surface  (^  =  wt)  and  incorporation  of 


t  t  t  t 


(a)  <b) 

Figure  14:  (a)  Structural  acoustic  system;  (b)  Piezoceramic  patches  driven  out-of-phase  to  create  pure 
bending  moments. 
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pressure  as  a  force  on  the  beam  then  yields  the  coupled  model 

Pf  d'2(t>< 


c2  dt2 


^dw  T 


V£e?u; 


+ 


+ 


[  PfVcf)- 
Jfi 

r  d2w  ,  r  „ 

Jr/SF^  +  kE-S*iSxi 

f  (9(t>  dWt\j 

Jr*  {at’1 -mV*1 


Td2w  d2rj 

EI-^—^-^—^dnf  + 


L 


d3w  62r\ 
CDl dx2dt  dx2^ 


(33) 


for  all  (£,  r/)  E  V.  The  space  of  test  functions  is  V  =  II1  (Q)  x  17q  (IT)  where  II1  (0)  denotes  the  quotient 
space  of  II [  (0)  over  constant  functions  and  17q  (To)  is  the  subset  II2 {To)  whose  elements  satisfy  the  fixed 
boundary  conditions  for  the  beam.  This  model  provides  a  baseline  for  comparison  since  it  incorporates 
the  full  physics  for  the  system. 

A  second  technique  used  to  model  the  acoustic  field  produced  by  the  vibrating  structure  is  to 
estimate  the  RMS  sound  power  radiated  from  the  structure  through  a  discretization  of  Raleigh’s  integral. 
Radiation  filters  are  then  constructed  to  incorporate  the  physical  mechanisms  modeling  the  structural 
acoustic  coupling.  These  filters  are  combined  with  a  state  variable  model  of  the  structure,  which  includes 
the  input  and  output  provided  by  piezoceramic  actuators  and  sensors,  to  provide  an  augmented  plant 
model.  Details  concerning  this  modeling  approach  for  structural  acoustic  systems  can  be  found  in  [19]. 
While  the  accuracy  of  the  approximations  to  the  Raleigh  integrals  degrades  near  boundaries,  the  method 
can  lead  to  significantly  lower-order  systems  than  obtained  through  a  full-order  approximation  of  (33). 
Furthermore,  the  prediction  of  sound  levels  solely  through  consideration  of  the  structural  vibrations 
permits  the  elimination  of  acoustic  sensors  and  promotes  the  use  of  smart  material  transducers  such  as 
piezoceramic  sensoriactuators.  Hence  the  smart  structure  capabilities  provided  by  these  materials  can 
be  utilized  in  structural  acoustic  applications. 


2.5  Control  Applications 

2.5.1  Infinite  and  Finite  Dimensional  Control  Problems 

The  problem  of  determining  optimal  voltages  to  patches  employed  as  actuators  in  structural  and  struc¬ 
tural  acoustic  systems  leads  to  the  formulation  of  infinite  dimensional  control  problems  with  the  models 
(25)  or  (33)  as  constraints.  Due  to  the  nature  of  the  patch  inputs  and  outputs,  these  control  problems 
typically  involve  unbounded  (discontinuous)  input  and  output  operators  and  are  quite  often  formulated 
as  boundary  control  problems.  The  analysis  of  the  linear  quadratic  regulator  (LQR)  problems  associated 
with  these  systems  and  the  determination  of  convergence  criteria  for  the  resulting  finite  dimensional 
systems  has  resulted  in  the  development  of  a  fairly  extensive  theory  for  linear  models.  Corresponding 
infinite  dimensional  theory  for  control  systems  with  nonlinear  or  hysteretic  inputs  is  far  less  complete. 
A  detailed  discussion  of  the  infinite  dimensional  LQR  theory  associated  with  smart  structures  is  beyond 
the  scope  of  this  paper  and  the  reader  is  referred  to  [3,  5,  9,  55]  for  details  concerning  the  state  of  this 
theory. 

In  subsequent  discussion,  we  focus  on  finite  dimensional  control  systems  of  the  form 

x(t)  =  Ax(t)  +  [B(u)](t)  +  g(t) 
y(t)  =  Cx(t ) 
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where  x  is  a  TV  x  1  state  vector,  y  denotes  observations  in  and  C  is  a  P  x  iV  observation  matrix 
whose  elements  are  specified  by  the  form  of  the  observations  (e.g.,  see  [9,  page  232]).  In  this  section,  we 
focus  on  the  linear  case  [B(u)](t)  —  Bu(t)  where  B  is  an  N  x  s  matrix  for  structures  with  s  actuator 
pairs.  We  then  consider  the  nonlinear  case  in  Section  3.5  in  the  context  of  magnetostrictive  actuators. 

2.5.2  Feedback  Control  Design 

We  consider  two  cases  for  the  exogenous  force  g(t ):  (i)  g(t)  =  0  which  can  be  used  when  controlling 
systems  subject  to  stochastic  or  modeling  uncertainties,  and  (ii)  _<y ( 0 )  =  g{r)  which  arises  when  modeling 
periodic  exogenous  inputs.  For  the  first  case,  the  LQR  problem  concerns  the  determination  of  a  control 
u  which  minimizes 


{(Qx(t) ,  x(t))  +  (Ru(t),u(t))}  dt 


(34) 


subject  to 


x(t)  —  Ax(t)  +  Bu(t) 
y(t)  -  Cx(t ) . 


(35) 


The  matrix  Q  is  chosen  to  satisfy  kinetic  and  potential  energy  criteria  and  R  weights  the  voltage 
to  the  s  patch  pairs  [4],  While  theoretically  attractive,  the  minimization  of  (34)  requires  full  state 
information  which  is  typically  not  available  in  applications.  This  necessitates  the  consideration  of  a 
dynamic  compensator  which  estimates  or  reconstructs  states  before  computing  a  feedback  control. 

As  detailed  in  [8,  9],  the  optimal  full-order  control  is  specified  by 


u(t)  —  —Kxc(t) 


(36) 


where  the  estimated  state  satisfies 

xc{t)  —  Axc(t)  +  Bu(t )  +  F  [y(t)  -  Cxc(t)\ 
xc{0)  =  xC0  . 

Here  K  and  F  denote  the  compensator  and  feedback  gains,  respectively.  We  note  that  K  and  F  are 
chosen  so  that  the  reconstruction  error  \y(t)  —  yc{t)  |  — >  0  as  t  — >  oo.  Under  usual  observability  and 
controllability  hypotheses  (see  [52]),  the  optimal  feedback  and  compensator  gains  are  given  by 

k  =  rt1bt  n 

F  =  PC'1!!-1 


(38) 


where  n  and  P  are  unique  nonnegative-definite  solutions  to  the  feedback  (regulator)  and  compensator 
(observer)  algebraic  Riccati  equations 

n  a  +  at  n  -  ubr-lbtu  +  q  =  o 

(39) 

PAt  +  AP-  PCTR~lCP  +  Q  =  0  , 

respectively.  As  was  the  case  with  the  matrices  Q  and  R.  the  matrices  Q  and  R  are  design  criteria  for 
the  specific  control  application  under  consideration. 
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For  case  (ii)  in  which  g  is  assumed  periodic  and  r  is  commensurate  with  all  frequencies  present  in 
the  exogenous  signal,  the  optimal  control  is  determined  by  minimizing 


Jt{u)  =  \  jo  {(Qz(f),z(f))  +  {Ru{t),u{t))} 


dt 


(40) 


subject  to 

x(t)  =  Ax(t)  +  Bu(t)  +  g(t)  ,  x(0)  =  x{t) 

y{t)  =  Cx{t )  . 

The  optimal  control  in  this  case  is  given  by 

u(t)  =  —kxc(t)  +  R~lBT z(t) 


(41) 


(42) 


where  the  reconstructed  state  is  specified  by  the  differential  equation 

xc(t)  =  Axc(t)  +  Bu(t)  +  F  \y(t)  -  Cxc{t)\  +  g(t ) 
£c(0)  =  xc(r) . 


and  the  tracking  variables  are  solutions  to  the  adjoint  equation 

z(t)  =  -[A  -  BK]Tz{t)  +  n g(t) 
z( 0)  =  z(t)  . 


(44) 


The  feedback  and  observer  gains  K  and  F  are  defined  in  (38)  and  II  is  the  solution  to  the  feedback 
Riccati  equation  in  (39). 

We  refer  the  reader  to  [8,  9]  for  details  concerning  the  original  infinite  dimensional  control  problem, 
criteria  which  guarantee  the  convergence  of  finite  dimensional  gains,  and  examples  illustrating  the  ex¬ 
perimental  implementation  of  the  control  method  for  a  circular  plate  with  surface-mounted  piezoceramic 
actuators. 


2.5.3  Adaptive  Feedforward  Control 

Currently  employed  adaptive  feedforward  algorithms  are  based  on  the  representation  of  the  plant  and 
control  as  finite  impulse  response  (FIR)  filters  whose  coefficients  are  determined  using  LMS  algorithms. 
This  discussion  provides  only  an  overview  of  the  methods  and  readers  are  referred  to  [19]  for  details 
concerning  feedforward  control  design  for  smart  or  adaptive  structures,  [62]  for  additional  analysis 
concerning  acoustic  systems,  and  [50,  89]  for  general  theory  of  adaptive  filters  and  feedforward  control 
design.  As  will  be  noted  in  this  summary,  feedforward  control  methods  are  based  on  superposition 
principles  and  hence  are  linear  by  nature.  Thus  they  are  restricted  to  linear  systems  or  systems  which 
have  been  suitably  linearized. 

Feedforward  control  design  for  structural  or  structural  acoustic  systems  which  employ  piezoceramic 
sensors  or  actuators  can  be  considered  in  two  steps:  (1)  Represent  the  plant  as  an  FIR  filter  and  (2) 
Construct  an  adaptive  filter  to  determine  the  control.  For  a  discrete  signal  u(k),  a  digital  filter  of  order 
L  +  1  for  the  plant  is 


L 

y{k)  =  E  h((k)u(k) 

e=o 


(45) 


23 


Figure  15:  System  identification  of  a  smart  structure  using  a  two  coefficient  FIR  filter  (after  [19]). 

where  h(k)  =  [ho (&),•••  ,  denotes  weights  which  must  be  determined  through  some  optimization 

procedure  (e.g.,  steepest  descent,  Newton’s  method,  LMS  algorithms).  As  detailed  in  [19],  the  minimiza¬ 
tion  of  a  quadratic  performance  functional  while  employing  instantaneous  estimates  of  the  correlation 
matrix  and  cross-correlation  vector  yields  the  LMS  algorithm 

h(k  +  1)  =  h(k)  +  nu(k)e(k)  (46) 

where  e(k)  =  d(k)  —  y(k)  and  //  is  chosen  so  that  filter  weights  approach  the  optimal  Wiener  filter  (see 
also  [88]).  A  schematic  of  the  system  identification  process  is  depicted  in  Figure  15.  The  transversal 
filter  (45)  is  used  both  to  represent  the  response  of  the  system  to  an  applied  control  and  to  filter  reference 
signals  applied  to  the  plant. 

The  second  step  is  the  determination  of  an  adaptive  FIR  filter  for  the  control.  Considering  first  the 
case  of  a  single  input,  the  control  is  represented  by  the  FIR  filter 

N 

u{k)  =  '^2wnr{k  —  n)  (47) 

n= 0 

where  r(k  —  n)  denotes  the  reference  signal  and  wn  are  filter  coefficients.  From  (45),  it  follows  that  the 
output  of  the  plant  in  response  to  the  control  input  is 

L  N 

Vu{k)  =  he  XI  wnr(k  ~n~  l). 

1=0  n= 0 

The  total  plant  output  is  then  represented  by  the  convolution 

y(k)  =  yr[k)-yu{k) 

N 

=  Vr[k)  -^2wnru{k  -n) 
n—1 
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Generalized  Plant 


Figure  16:  Filtered-x  feedforward  control  algorithm  as  applied  to  the  smart  material  structure  (af¬ 
ter  [19]). 

where 

L 

ru(k  —  n)  =  ^  htr(k  —  n  —  i) . 
t=  o 

The  signal  ru(k  —  n)  is  known  as  the  filtered  reference  signal  and  from  (45),  it  is  observed  to  be  the 
output  of  the  plant  in  response  to  the  reference  signal  r(k).  For  implementation,  this  signal  can  be 
computed  off-line  using  the  LMS  algorithm  (46)  to  determine  the  plant  coefficients  he .  This  is  the 
filtered-x  step  in  the  control  computation. 

The  control  coefficients  are  determined  by  minimizing  the  cost  functional 

J  =  E[y2(k )] 

using  the  approximation  techniques  of  Widrow  and  Stearns  [89].  This  yields  the  adaptive  filtered-x 
algorithm 

w(k  +  1)  =  w(k)  +  p,ru(k)y(k) 

where  /}  is  a  second  gain  parameter  which  regulates  the  stability  and  convergence  of  the  controller.  The 
structure  of  the  algorithm  is  depicted  in  Figure  16. 

As  detailed  in  [19],  the  multi-input  multi-output  case  is  analyzed  similarly  and  yields  the  MIMO 
filtered-x  LMS  algorithm  with  P  actuators 

w(&  +  1)  =  w  (k)  +  firu{k)y{k) 

where 

W (fc)  =  [wj (k),W2  (ft),*w  ,Wp{k)] 
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and 

wp{k)  =  [wp,i{k),wp,2{k),- ■  ■  ,wp,N{k)]T  ,  p  =  1,2,  -  ■  ■  ,  P . 

Details  regarding  the  construction  of  the  matrix  ru(k)  for  the  filtered  reference  signal  can  be  found 
in  [19]. 

2.5.4  Hybrid  Control 

Hybrid  control  methods  can  be  designed  to  combine  the  strengths  of  feedback  and  feedforward  formula¬ 
tions  to  improve  controller  performance  in  systems  subjected  to  both  stochastic  and  harmonic  exogenous 
forces  as  well  as  systems  containing  modeling  uncertainties.  Furthermore,  the  shifting  of  plant  poles 
through  the  incorporation  of  a  feedback  compensator  in  the  control  path  can  significantly  improve 
convergence  times  for  the  feedforward  component  of  the  system.  Details  and  software  illustrating  the 
implementation  of  such  hybrid  methods  can  be  found  in  [19]. 

To  indicate  the  performance  capabilities  of  such  hybrid  approaches,  we  consider  the  beam  illustrated 
in  Figure  13  with  a  pair  of  piezoceramic  actuators.  A  600  Hz  exogenous  force  p  was  applied  uniformly 
along  the  beam  to  simulate  a  periodic  driving  pressure  field.  The  beam  was  assumed  to  have  length 
l  =  .46  m,  width  b  =  .02  m  and  thickness  h  =  .0016  m,  and  the  patches  extended  from  x\  =  .2  m  to 
X2  =  .25  m.  The  density,  stiffness  and  damping  parameters  in  (22)  and  (23)  were  specified  to  be  p(x)  = 
.093 +.433%^  (a;)  kg/m,  EI(x)  =  .491  +  .  793;^^  (a;)  Nm2  and  cdI(x)  =  .649  x  10-5 +  1.255  x  10 ~5Xpzt(x) 
Nm2sec  in  accordance  with  physical  parameters  identified  for  similar  structures. 

A  hybrid  controller  comprised  of  the  feedback  law  (36)  implemented  in  parallel  with  the  feedforward 
law  (47)  was  used  to  determine  controlling  voltages  to  the  patches  given  point  measurements  of  the 
velocity  at  the  point  x  =  ‘M/5.  The  controlled  velocity  and  adaptation  of  the  weights  in  the  feedforward 
filter  are  plotted  in  Figure  17.  It  is  observed  that  by  0.1  seconds  the  weights  have  converged  and  by 
0.2  seconds,  vibrations  have  been  effectively  attenuated.  When  considered  in  combination  with  cited 
references,  this  illustrates  aspects  of  representative  linear  control  systems  which  are  currently  employed 
for  smart  structures  which  utilize  piezoceramic  actuators  in  their  linear  range. 


(a)  (b) 

Figure  17:  (a)  Beam  velocity  at  the  point  x  =  M/5,  (b)  Adaptation  of  the  filter  weights  to  accommodate 
the  harmonic  exogenous  force  to  the  beam. 
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3  Magnetostrictive  Materials 


Magnetostrictive  transducers  are  increasingly  employed  in  automotive  and  industrial  applications  which 
require  the  generation  of  massive  forces  or  strains,  or  the  sensing  of  magnetic  or  elastic  quantities  such 
as  torques  or  strains.  These  transducer  materials  rely  on  the  magnetostrictive  effect  which  is  inherent 
to  all  ferromagnetic  materials  but  is  sufficiently  large  in  compounds  such  as  nickel  or  Terfenol-D  to 
warrant  their  use  as  high  performance  actuators  and  sensors.  For  such  compounds,  the  application  of 
a  magnetic  field  produces  strains  in  the  material  while  conversely,  strains  produced  by  external  forces 
generate  changes  in  the  magnetic  state.  Because  the  magnetostriction  is  an  inherent  property  of  the 
materials,  their  performance  does  not  degrade  over  time  as  can  be  the  case  with  piezoelectric  compounds. 

To  motivate  modeling  and  control  issues,  we  first  outline  certain  applications  which  utilize  magne¬ 
tostrictive  actuators  and  sensors.  We  then  summarize  pertinent  material  characteristics  as  a  prelude 
to  the  development  of  constitutive  and  hysteresis  models  and  their  incorporation  in  system  models.  As 
in  the  previous  section  on  piezoelectric  and  electrostrictive  transducers,  we  conclude  with  prototypical 
control  designs  and  directions  for  future  research. 

3.1  Magnetostrictive  Transducer  Applications 

The  effects  of  magnetostriction  have  been  recognized  since  the  1840’s  and  devices  utilizing  these  prop¬ 
erties  have  been  in  existence  since  the  1860’s  [38].  In  the  first  half  of  this  century,  applications  uti¬ 
lizing  nickel  and  other  magnetostrictive  compounds  included  telephone  receivers,  torque  sensors,  and 
hydrophones  [24],  In  1978,  A.E.  Clark  and  co-workers  developed  Terfenol-D  (Tb^Dyi-^Fe)  which  is 
classified  as  a  “giant”  magnetostrictive  alloy  due  to  its  capability  for  generating  strains  in  excess  of 
1000  n L/L.  Applications  which  utilize  these  alloys  include  active  vibration  and  noise  control,  microp¬ 
ositioning  in  high  force  regimes,  medical  and  industrial  ultrasonics,  non-contact  torque  sensors,  and 
tuned  vibration  absorbers.  The  reader  is  referred  to  [24,  69]  for  a  comprehensive  discussion  of  these 
applications. 

The  design  of  magnetostrictive  transducers  is  rapidly  evolving  as  the  capabilities  of  these  materials 
become  better  understood.  A  common  design  for  numerous  structural  applications,  which  motivates 
pertinent  modeling  and  control  issues,  is  the  piston  actuator  depicted  in  Figure  18.  Strains  in  the 
Terfenol-D  rod  are  produced  when  fields  generated  by  current  through  the  solenoid  cause  the  rotation 
of  magnetic  moments  in  the  manner  illustrated  in  Figure  19.  The  sensor  capabilities  of  the  transducer  are 
produced  by  the  opposite  effect  in  which  stresses  to  the  rod  realign  moments  which  produces  subsequent 
changes  in  the  magnetization.  The  prestress  mechanism  serves  two  roles:  it  increases  the  percentage 
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Figure  18:  Cross  section  of  a  prototypical  Terfenol-D  magnetostrictive  transducer. 
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(c)  e  =  A. s 


Figure  19:  Magnetic  domains  in  the  Terfenol-D  rod;  (a)  Orientation  of  unstressed  rod  in  absence  of 
applied  magnetic  field,  (b)  Orientation  of  pre-stressed  rod  with  no  applied  field,  (c)  Orientation  of 
pre-stressed  rod  when  field  is  applied  in  direction  of  longitudinal  rod  axis. 

of  moments  which  are  oriented  perpendicular  to  the  rod  axis  and  it  maintains  the  rod  in  a  state  of 
compression.  The  bias  fields  required  to  attain  bidirectional  strains  and  improved  linearity  are  provided 
either  by  a  surrounding  permanent  magnet  or  through  a  prescribed  DC  current  to  the  solenoid.  In  spite 
of  the  fact  that  the  permanent  magnet  adds  bulk,  it  is  typically  employed  since  it  reduces  ohmic  heating 
to  the  coil  and  provides  additional  flexibility  for  shaping  the  flux  path  to  optimize  performance. 

The  manner  through  which  the  transducer  can  be  employed  in  acoustic  and  structural  systems 
is  illustrated  in  Figure  20.  As  indicated  in  (a)  and  (b),  the  transducers  can  be  used  to  both  sense 
vibrations  in  structures  or  industrial  equipment  and  provide  the  force  or  strain  inputs  required  to 
attenuate  vibrations.  A  current  area  of  development  is  the  utilization  of  magnetostrictive  materials  for 
non-contact  torque  sensing.  Figure  20(c)  illustrates  the  coupling  of  the  transducer  with  an  ultrasonic 
horn  for  applications  ranging  from  the  cleaning  of  intricate  or  inaccessible  machinery  to  the  catalysis  of 
chemical  reactions  [24], 

A  second  class  of  magnetostrictive  transducers  includes  flextensional  devices  analogous  to  those 
described  in  Section  2.1,  inchworm  devices  and  rotational  motors.  The  flextensional  designs  are  similar 


(c) 


Figure  20:  Applications  utilizing  magnetostrictive  transducers,  (a)  Vibration  sensing  and  attenuation, 
(b)  Torque  sensing  and  vibration  attenuation  in  a  milling  machine,  (c)  Ultrasonic  horn  (after  [24]). 
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(b) 


Figure  21:  Hybrid  magnetostrictive-piezoceramic  transducers,  (a)  Inchworm  actuator  (after  [59])  and 
(b)  Rotational  motor  (after  [86]). 


to  those  depicted  in  Figure  4  but  utilize  a  magnetostrictive  core  rather  than  PZT  or  PMN  stacks. 
The  magnitude  of  generated  forces  in  combination  with  the  dual  actuating  and  sensing  capabilities 
of  the  materials  makes  the  transducers  potentially  advantageous  in  hydrostatic  applications  including 
sonar  transduction  and  underwater  imaging.  Hybrid  devices  employing  magnetostrictive  drive  rods  and 
piezoceramic  clamps  are  depicted  in  Figure  21.  In  addition  to  the  forces  provided  by  the  magnetostrictive 
materials,  the  intrinsic  nature  of  the  magnetomechanical  coupling  provides  the  possibility  for  eliminating 
the  variability  and  aging  associated  with  conventional  hydraulic  and  electromagnetic  motors. 

A  comparison  of  magnetostrictive  transducers  with  the  piezoelectric  and  electrostrictive  devices 
discussed  in  Section  2.1  indicates  that  while  they  exhibit  superior  performance  in  certain  applications, 
they  are  bulkier  than  their  ferroelectric  counterparts  due  to  the  required  solenoid  and  housing.  Hence 
they  are  currently  employed  in  regimes  where  giant  forces  or  strains  are  required  but  weight  is  not  an 
issue.  The  range  of  applications  will  almost  surely  increase,  however,  as  materials  such  as  magneto¬ 
composites  are  perfected. 

While  there  are  variations  among  magnetostrictive  transducer  designs,  certain  modeling  and  control 
issues  are  common  to  all  applications.  Hysteresis  is  present  to  some  degree  in  all  magnetostrictive  ma¬ 
terials  and  the  constitutive  relations  between  the  magnetization  and  strain  exhibits  quadratic  behavior 
at  low  drive  levels  and  saturates  at  high  drive  levels.  Hence  biasing  fields  or  stresses  are  often  applied 
to  produce  bidirectional  strains  and  improve  linearity.  In  addition  to  quantifying  the  nonlinear  and 
hysteretic  relations  between  the  field,  magnetization  and  strain,  it  is  necessary  to  characterize  the  elec¬ 
tromagnetic  relations  between  the  current  input  to  the  solenoid  and  the  field  generated  at  various  points 
in  the  rod,  the  frequency-dependent  losses,  the  coupled  magnetomechanical  dynamics  of  the  structure, 
and  thermal  effects.  While  aspects  of  this  have  been  accomplished  for  specific  transducer  configurations, 
comprehensive  models  which  incorporate  all  of  these  mechanisms  are  lacking.  In  Sections  3.3  and  3.4, 
we  summarize  certain  constitutive  and  hysteresis  models  and  illustrate  their  incorporation  in  system 
models  for  smart  structures.  The  development  of  commensurate  control  methods  is  then  outlined  in 
Section  3.5. 

3.2  Material  Behavior 

The  hysteresis  and  nonlinear  constitutive  behavior  of  magnetostrictive  materials  are  due  in  part  to  the 
domain  structure  inherent  to  ferromagnetic  compounds.  We  summarize  here  material  attributes  which 
directly  pertain  to  the  performance  of  the  compounds  in  transducer  design  both  to  illustrate  intrinsic 
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properties  of  the  materials  and  to  motivate  the  models  summarized  in  Section  3.3.  We  note  that 
domain  mechanisms  in  ferromagnetic  materials  are  analogous  in  many  ways  to  the  domain  attributes 
which  produce  hysteresis  in  ferroelectric  materials  (Section  2.3)  and  when  possible,  we  try  to  consider 
models  and  control  designs  which  encompass  both  sets  of  compounds. 

The  description  of  crystalline  attributes  of  the  materials  is  facilitated  by  vector  conventions  in  which 
planes  are  denoted  by  round  brackets  and  directional  indices  are  represented  by  square  brackets.  Hence 
the  faces  of  a  cube  are  denoted  by  (100),  (010),  (001),  (100),  (010),  (001)  while  the  vertices  are  specified 
by  [100],  [010],  [001],  [100],  [010],  [001].  In  both  cases,  1,  indicates  a  negative  direction.  Finally,  pointed 
brackets  are  used  to  summarize  an  entire  set  of  indices  (e.g.,  (100)  represents  the  vertices  of  the  cube). 

3.2.1  Ferromagnetic  Domain  Structure  and  Hysteresis 

For  current  smart  structure  applications,  the  most  widely  employed  magnetostrictive  material  is  Terfenol- 
D.  As  detailed  in  [16]  and  depicted  in  Figure  22,  Terfenol-D  crystals  are  grown  in  dendrite  sheets  oriented 
in  the  [112]  direction.  At  room  temperature,  the  easy  axes  for  Terfenol-D  lie  approximately  in  the  (111) 
direction  and  the  greatest  strains  occur  when  the  magnetization  M  rotates  from  [111]  to  [111]. 

The  changes  in  magnetism  which  result  from  an  applied  field  H  are  primarily  due  to  two  mechanisms: 
the  rotation  of  moments  and  the  movement  of  domain  walls.  To  illustrate,  consider  a  demagnetized 
compound  which  is  subjected  to  a  magnetic  field  oriented  in  the  [112]  direction  as  shown  in  Figure  23. 
At  low  field  levels,  the  primary  change  in  magnetization  is  due  to  domain  wall  motion  so  that  favorably 
oriented  domains  are  enlarged.  As  the  field  is  increased,  domains  rotate  to  orient  with  the  easy  [111] 
axis.  This  produces  a  burst  region  in  the  H-M  or  H-e  curve  in  which  small  changes  in  field  produce 
large  changes  in  magnetization  or  strain  (see  Figure  24).  From  a  design  perspective,  biasing  into  this 
region  can  be  used  to  optimize  performance.  In  the  final  stage  depicted  in  Figure  23,  the  material  acts 
as  a  single  domain  as  dipoles  rotate  from  the  easy  axis  into  the  direction  of  the  applied  field.  This 
produces  the  saturation  observed  in  Figure  24. 

For  low  field  deviations  about  a  bias  field  level,  the  relation  between  the  field  and  magnetization 
is  approximately  linear  as  illustrated  in  Figure  24a.  At  higher  field  levels,  however,  magnetostrictive 
materials  typically  exhibit  hysteresis  of  the  form  depicted  in  Figure  24b.  As  with  ferroelectric  compounds 


Figure  22:  Orientation  of  Terfenol-D  crystals. 


30 


Figure  23:  Magnetization  process  due  to  an  applied  field  H.  (a)  Demagnetized  state  and  (b)  growth 
of  domains  due  to  domain  wall  motion,  (c)  Rotation  of  domains  to  the  easy  axis  and  (d)  Rotation  of 
domains  to  align  with  the  applied  field. 


(see  Section  2.2),  this  hysteresis  can  be  due  to  the  irreversible  impediment  of  domain  wall  movement  by 
inclusions  or  pinning  sites  in  the  material.  Details  regarding  the  sources  of  hysteresis  in  ferromagnetic 
materials  can  be  found  in  [10,  22,  44], 

Strains,  or  magnetostriction,  are  generated  by  the  material  when  dipoles  rotate  to  align  with  an 
applied  field.  For  general  configurations,  the  magnetomechanical  coupling  which  produces  these  strains 
is  highly  complex  and  dependent  upon  factors  such  as  the  applied  stress  and  crystalline  anisotropies. 
However,  in  the  case  of  materials  in  which  a  prestress  perpendicular  to  the  dipole  direction  is  sufficiently 


(a)  (b) 

Figure  24:  (a)  Anhysteretic  magnetization  curve,  (b)  Ferromagnetic  hysteresis  present  at  high  drive 
levels. 
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large  to  dominate  crystalline  anisotropies,  the  strains  are  due  primarily  to  domain  rotation  in  the  manner 
depicted  in  Figure  19.  For  this  regime,  the  magnetostriction  A  can  be  modeled  by  the  quadratic  relation 


(48) 


where  As  and  Ms  respectively  denote  the  saturation  magnetostriction  and  magnetization.  As  discussed 
in  [13,  44],  the  prestress  levels  which  optimize  transducer  performance  are  often  of  a  magnitude  so  that 
stress  anisotropies  dominate  crystalline  anisotropies  and  the  relation  (48)  adequately  models  the  free 
strains  generated  by  the  material. 


3.3  Transducer  Models 

A  large  number  of  modeling  techniques  have  been  employed  to  characterize  the  behavior  of  magne- 
tostrictive  materials  in  various  operating  regimes.  These  include  linear  and  higher-order  anhysteretic 
models  as  well  as  micromagnetic,  Preisach,  domain  rotation  and  domain  wall  models  for  the  hysteresis 
present  at  higher  drive  levels.  We  outline  here  models  which  have  been  employed  for  control  design 
and  direct  the  reader  to  various  references  for  more  comprehensive  perspectives  concerning  individual 
modeling  techniques. 


3.3.1  Anhysteretic  Models 

For  low  drive  level  applications  in  which  hysteresis  is  minimal,  anhysteretic  (hysteresis-free)  models 
analogous  to  those  employed  for  piezoceramic  and  electrostrictive  compounds  have  been  adopted  for 
magnetostrictive  materials.  In  most  control  applications,  a  bias  field  produced  either  by  a  surrounding 
permanent  magnet  (see  Figure  18)  or  DC  field  to  the  solenoid  is  used  to  attain  bidirectional  strains  and 
improve  linearity.  For  low  field  levels,  the  linear  model 

B-i  =  dijk&jk  +  HijHj 

H  / 

dij  ^ijki^ki  +  dhijHfc 

is  often  employed.  A  comparison  with  (1)  reveals  the  similarity  with  the  linear  piezoelectric  model. 
The  stress,  strain  and  compliance  tensors  cr,  e  and  sH  are  identical  to  those  employed  in  Section  2.3 
while  the  tensor  d  is  analogous  to  the  piezoelectric  strain  tensor.  Finally,  jp  denotes  the  permeability 
at  constant  stress.  The  model  (49)  has  the  advantage  of  simplicity  and  linearity  but  does  not  incorpo¬ 
rate  the  nonlinear  constitutive  behavior  or  hysteresis  present  in  many  high  drive  level  smart  structure 
applications. 

Saturation  phenomena,  temperature  effects,  and  the  quadratic  relation  between  the  magnetization 
and  free  strains  summarized  in  (48)  are  incorporated  in  the  higher-order  anhysteretic  model 

Hk  =  -2QjjktMtatj  +  — arctanh  (\M\  jMs )  +  P%  AT 

ki\M\  (50) 

eij  —  sijHake  +  QijkfMkMf,  +  Oifj  AT 

of  Duenas,  Hsu  and  Carman  [28].  Here  QT  is  a  fourth  rank  tensor  of  magnetostrictive  coefficients,  T  is 
the  temperature,  aM  is  a  coefficient  of  thermal  expansion,  Pa  is  a  pyromagnetic  coefficient  and  k\  is  a 
constant.  In  each  case,  the  superscript  denotes  quantities  which  are  fixed.  A  comparison  with  (4)  for  T 
given  by  (9)  illustrates  that  (50)  is  analogous  to  the  quadratic  electrostriction  model  in  which  the  Ising 
spin  relation  (7)  is  used  to  characterize  the  anhysteretic  polarization. 
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The  higher  order  model  quantifies  certain  aspects  of  the  nonlinear  constitutive  behavior  but  can 
incorporate  hysteresis  only  through  the  use  of  complex  model  parameters.  Techniques  for  modeling  the 
hysteresis  in  the  magnetization  77  in  the  first  (converse)  relation  in  (50)  are  addressed  next. 

3.3.2  Hysteresis  Models 

Hysteresis  models  for  magnetostrictive  materials  range  from  micromechanical  models  to  phenomenologi¬ 
cal  characterizations.  Micromechanical  models  quantify  the  coupled  changes  in  magnetization  and  strain 
through  first  principles  formulations  of  the  elastic  and  magnetostatic  states  [11,  42],  While  this  can 
produce  highly  accurate  material  models,  the  large  number  of  required  parameters  currently  precludes 
the  real-time  implementation  of  controllers  which  employ  these  models.  Preisach  models  for  magnetic 
and  magnetostrictive  materials  have  been  developed  by  numerous  researchers  [1,  70,  76],  and  the  advan¬ 
tages  and  disadvantages  of  this  approach  are  similar  to  those  discussed  for  ferroelectric  Preisach  models 
(see  Section  2.3).  Domain  and  domain  wall  models  include  classical  and  modified  anisotropy  domain 
rotation  models  [46,  84]  and  the  domain  wall  theory  of  Jiles  and  Atherton  [45]  which  motivated  the 
ferroelectric  hysteresis  model  discussed  in  Section  2.3  for  piezoceramic  and  electrostrictive  models. 

As  detailed  in  Section  2.3,  the  domain  wall  theory  models  hysteresis  by  quantifying  the  energy 
required  to  bend  and  translate  domain  walls  pinned  at  inclusions  inherent  to  the  material.  This  yields 
irreversible  and  reversible  components  whose  sum  is  the  total  polarization  in  ferroelectric  materials  or 
magnetization  in  ferromagnetic  materials. 

From  the  theory  of  Jiles  and  Atherton  [45],  the  irreversible  magnetization 

.  dMirr 

Mirr  =  Man  -  k 

dHe 

is  specified  by  quantifying  losses  from  the  minimum  energy  anhysteretic  magnetization  Man.  The 
effective  field  is  given  by  He  —  H  +  aM  and  the  anhysteretic  magnetization  is  modeled  by  the  Langevin 
equation 

M0„  =  Ms[coth(f)-(£)[ 

where  a  and  a  are  parameters  specifying  the  slope  of  the  anhysteretic  curve.  The  parameter  k  quan¬ 
tifies  the  average  energy  required  to  change  moment  orientations  and  hence  break  pinning  sites.  The 
formulation  of  the  irreversible  magnetization  in  terms  of  the  applied  field  yields 

d,Mvrr  _ ~  Man  Mvrr 

dH  kS  —  a  ( Man  —  Mirr ) 

where  S  —  sign (17)  and  S  is  defined  in  a  manner  analogous  to  (12). 

The  Jiles-Atherton  theory  characterizes  the  magnetization  due  to  the  bending  of  domain  walls 
through  the  algebraic  relation 

Mrev  —  c(Man  Mirr) 

where  c  quantifies  reversible  domain  wall  effects.  The  total  magnetization  is  then  specified  by 

77  —  77^rr  T  77rei, .  ( 5 1 ) 

We  note  that  the  formulation  of  this  model  is  identical  to  the  ferroelectric  hysteresis  model  (14).  This 
permits  both  a  unified  theory  for  characterizing  hysteresis  in  the  two  classes  of  materials  and  the 
development  of  control  algorithms  which  accommodate  both  classes  of  transducers. 
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An  alternative  formulation  of  the  model  which  facilitates  inversion  can  be  obtained  by  consolidating 
the  irreversible  and  reversible  magnetizations  to  obtain  the  differential  equation 


dM 

~dH 


E{H,  M) 


M(Hq)  =  M0 


(52) 


where 


F{H,  M) 


1 


1  +  cMsa§§(^) 


6  [Ms£  (^±f^)  -  M]  dC  (H  +  aM 

kd-a  [1 MSC  ( toM)  -  M]  °Ms  dH  \  a 


The  Langevin  function  is  defined  by 

C(z)  =  coth(z) - . 

The  time-dependent  dynamics  of  the  magnetization  can  then  be  specified  through  the  chain  rule 


dM 

dt 


E{H,M) 


dH 

dt 


(53) 


We  note  that  in  the  formulation  outlined  here,  frequency-dependent  eddy  current  losses  are  ignored  so 
the  relation  (53)  should  be  employed  in  regimes  where  these  effects  are  minimal. 

The  magnetization  model  (51)  or  (52)  can  be  combined  with  the  constitutive  relation  (50)  to  quantify 
the  strains  generated  by  the  transducer  when  hysteresis  is  significant.  Experimental  validation  of  this 
model  for  a  prototypical  transducer  of  the  type  depicted  in  Figure  18  can  be  found  in  [25,  26]. 

The  hysteresis  model  in  the  form  (52)  is  also  amenable  to  inversion  through  consideration  of  the 
complementary  differential  equation 


dM-1  _  1 

dH  ~  Y{aF\h) 

M-\H0)  =  Mq-1 


The  experimental  validation  of  the  inversion  procedure  is  detailed  in  [83]  and  numerical  examples 
describing  its  role  in  control  design  will  be  outlined  in  Section  3.5. 


3.4  System  Models 

To  illustrate  the  effect  of  the  nonlinear  input  on  a  smart  structure  which  employs  magnetostrictive 
transducers,  we  consider  a  prototypical  beam  with  end-mounted  transducers  as  depicted  in  Figure  20a. 
The  description  of  the  beam  was  provided  in  Section  2.4  and  we  focus  here  solely  on  the  Terfenol-D 
transducers  which  are  driven  diametrically  out-of-phase  to  produce  pure  bending  moments  in  the  beam. 
As  detailed  in  [14],  this  setup  has  been  employed  in  experiments  to  ascertain  capabilities  and  properties 
of  the  transducers. 

For  modeling  purposes,  the  cross-sectional  area  and  Young’s  modulus  for  the  Terfenol-D  rod  are 
denoted  by  Amag  and  EH  while  the  length  of  the  connecting  rod  is  tr.  As  detailed  in  [75],  the  transverse 
displacements  of  the  beam  are  modeled  by  (21)  with  the  density  p.  stiffness  El  and  damping  coefficient 
cpl  modified  to  incorporate  material  properties  of  the  Terfenol  rod  and  connecting  bar.  For  a  system 
in  which  the  magnetization  is  biased  to  obtain  bidirectional  strains,  the  quadratic  relation  (48)  yields 
the  free  strains 

3A 

e(t)  =  ^[M2(t)+2M(t)Ms\. 
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The  external  moment  generated  by  the  transducers  is  then 

Mext(t,  x)  =  JCM[M2(t)  +  2 M{t)Ms]Xrod{x)  (55) 

where  JCM  =  ( 3 AiS /M2 ) Amag E 11  (h/2  +  ir)2  and  the  magnetization  M  is  modeled  by  the  hysteresis 
relation  (51)  or  (52). 

The  formulation  of  the  model  (21)  in  weak  form  and  approximation  using  the  Galerkin  expansion 
(28)  yields  the  first  order  system 


y(t)  =  Ay(t)  +  [£(«)](*)  +  g{t) 

y{  o)  =  yo 

where  the  nonlinear  input  term  has  the  form 

[B{u)](t)  =  [. M2(H{t ))  +  2 M(H{t))Ms\  b. 
The  system  matrices  and  vectors  y,  A,r  and  b  are  defined  in  (30). 


(56) 


3.5  Control  Formulation 

There  are  essentially  three  options  for  designing  control  methods  for  systems  which  employ  magne- 
tostrictive  transducers.  The  first  is  to  employ  the  linear  constitutive  relations  (49)  which  yields  a  linear 
system.  This  provides  a  setting  in  which  feedforward  and  feedback  methods  of  the  type  illustrated  in 
Section  2.5  can  be  employed,  but  is  accurate  only  at  low  drive  levels.  The  second  option  is  to  determine 
an  inverse  model  for  the  transducer  which  either  exactly  or  approximately  compensates  for  the  actuator 
hysteresis.  This  strategy  will  be  illustrated  here  for  magnetostrictive  transducers  and  in  the  context  of 
Preisach  models  for  shape  memory  alloys  in  Section  4.  The  final  option  is  to  design  a  nonlinear  control 
law  for  the  system.  This  is  done  in  [77]  for  the  structural  system  described  here  and  the  examples  in 
that  investigation  provide  a  baseline  for  determining  the  effectiveness  of  the  feedback  controls  which 
incorporate  an  inverse  compensator. 


3.5.1  Control  Example 

We  consider  two  techniques  for  formulating  a  linear  control  problem  for  the  structure  and  transducer 
setup  described  in  Section  3.4.  In  the  first  case,  we  simply  linearize  the  operator  B(u)  about  a  given  field 
level  u  =  H.  This  technique  is  commonly  employed  when  magnetostrictive  transducers  are  employed 
at  low  drive  levels  with  magnetic  biases  and  is  applicable  only  when  the  transducer  is  operated  in  the 
nearly  linear  range.  The  second  technique  is  based  upon  the  filtering  of  control  inputs  using  the  inverse 
compensator.  This  control  method  is  applicable  throughout  the  operational  range  of  the  transducer. 
Details  regarding  this  example  are  provided  in  [78]. 

The  spatially  uniform  exogenous  force 

g(t,x)  =  go  [sin(147rf)  +  sin(267rf)] 

was  applied  throughout  the  time  interval  [0, 2.5]  to  simulate  a  periodic  pressure  field  on  the  beam.  The 
uncontrolled  trajectories  at  the  point  x  =  3£/5  for  g0  =  1  are  plotted  in  Figures  25b  and  26a  while 
the  trajectories  obtained  with  g0  =  100  are  plotted  in  Figures  25d  and  26b  Both  cases  exhibit  a  beat 
phenomenon  due  to  the  close  proximity  of  the  7  Hz  driving  frequency  with  the  6.1  Hz  natural  frequency 
of  the  beam. 
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Figure  25:  Feedback  of  the  linear  law  (42)  into  the  nonlinear  system  (56).  Relationship  between  magnetic 

field  H  and  magnetization  M  for  (a)  go  =  1  and  (c)  go  =  100.  Uncontrolled  ( - )  and  controlled  (  ) 

beam  trajectories  at  x  =  3^/5  for  (b)  go  =  1  and  (d)  go  =  100. 


An  open  loop  control  input  was  computed  using  (42)  with  full  state  measurements  xc  =  x,  and 
applied  to  the  hysteretic  system  (56)  at  time  t  =  0.45.  The  resulting  controlled  beam  trajectories  at 
the  x  are  plotted  in  Figures  25b  and  25d  and  the  relationship  between  the  input  magnetic  field  H  =  u 
and  the  resulting  magnetization  M  for  the  two  cases  are  given  in  Figures  25a  and  c.  The  results  for 
go  =  1  illustrates  that  the  method  is  very  effective  at  low  drive  levels  where  the  linear  model  is  accurate. 
At  the  high  drive  levels  at  which  the  actuators  are  advantageous,  however,  the  input  B(u)  introduces 
significant  hysteresis  which  acts  as  a  time  delay  to  the  system.  The  result  is  a  loss  in  control  authority 
which  is  sufficiently  severe  to  produce  controlled  beam  trajectories  which  are  larger  in  magnitude  at 
certain  times  than  the  uncontrolled  trajectories. 

To  obtain  a  control  based  on  the  inverse  compensator,  we  then  employed  the  input  operator  B  =  b 
where  b  is  defined  in  (30).  This  yields  the  prescribed  input  u^.  To  obtain  the  optimal  control  for  the 
original  system  (56),  the  signal  is  then  filtered  through  the  inverse 


v(t)  =  M~l  (Ms  + 


(57) 


before  input  to  the  transducer.  An  overview  of  this  inversion  process  is  illustrated  in  Figure  12  in  the 
context  of  a  Preisach  model.  We  note  that  while  this  requires  the  on-line  integration  of  (54)  to  obtain 
M-1,  the  operation  is  significantly  less  expensive  than  the  online  integration  of  (44)  to  obtain  z(t)  or  the 
integration  of  (43)  required  to  obtain  state  estimates  xc(t)  when  full  state  information  is  unavailable. 


36 


(a)  (b) 

Figure  26:  Uncontrolled  beam  trajectory  ( - )  at  the  point  x  =  37/5  and  controlled  trajectory  (  ) 

obtained  with  the  filtered  input  (57);  (a)  go  =  1  and  (b)  go  =  100. 


The  controlled  trajectories  obtained  in  this  manner  for  go  =  1  and  go  =  100  are  respectively  plotted 
in  Figure  26.  In  this  case,  there  is  no  loss  of  control  authority  at  high  drive  levels,  thus  illustrating  that 
this  approach  provides  a  significant  advantage  over  the  control  design  based  on  the  linearized  system. 
Finally,  a  comparison  with  the  examples  in  [77]  illustrates  that  the  performance  of  the  controller  based 
on  inverse  compensation  is  comparable  to  that  obtained  when  the  functional  (40)  is  minimized  subject 
to  the  original  system  (56)  but  at  a  fraction  of  the  computational  cost. 

4  Shape  Memory  Alloys 

Shape  memory  alloys  (SMA)  comprise  a  class  of  metallic  compounds  which  convert  thermal  energy  into 
mechanical  energy.  While  the  underlying  mechanisms  are  complex,  the  generation  of  strains  in  response 
to  changes  in  temperature  or  stress  can  be  described  as  due  to  crystalline  transformations  from  a  high 
temperature  austenite  phase  to  a  low  temperature  martensite  phase  and  vice-versa.  For  the  most  com¬ 
monly  employed  material  Nitinol  (Ni:  nickel,  ti:  titanium,  nol:  Naval  Ordinance  Laboratory),  this  can 
lead  to  the  generation  of  recoverable  strains  on  the  order  of  6%.  Because  the  crystalline  transformations 
are  driven  by  changes  in  temperature,  the  original  shape  memory  alloys  had  a  low  bandwidth  which 
limited  their  utility  for  high  frequency  control  applications.  Their  potential  for  shape  modification  or 
control  is  significant,  however,  due  to  the  magnitude  of  the  generated  strains.  Furthermore,  the  devel¬ 
opment  of  hybrid  shape  memory  alloys  has  produced  transducer  materials  with  improved  bandwidth 
(e.g.,  see  [43]),  and  applications  utilizing  these  compounds  should  continue  to  evolve. 

The  nature  of  the  martensite  to  austenite  phase  transition  and  the  associated  change  in  shape  of 
the  SMA  produces  shape  memory  effects  (SME)  which  provide  the  materials  with  their  capabilities  as 
distributed  actuators.  The  phenomenological  behavior  of  the  phase  transition  is  illustrated  in  Figure  27a. 
The  cooling  of  the  SMA  from  the  high  temperature  austenite  phase  to  the  low  temperature  martensite 
phase  is  characterized  by  a  phase  transition  which  starts  at  a  temperature  Ms  and  is  finished  at  Mf. 
The  opposite  phase  transition  occurs  as  the  material  is  heated  through  the  temperatures  As  and  Aj. 

To  utilize  the  associated  shape  memory  effects,  the  SMA  are  stress-loaded  at  temperatures  T  <  As 
in  the  martensite  phase.  This  produces  the  residual  strain  er  depicted  in  Figure  27b.  These  residual 
strains  are  recovered  by  heating  the  compound  past  the  austenite  finish  temperature  Aj.  For  details 
concerning  this  process,  readers  are  referred  to  [9,  31,  56]. 
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Figure  27:  (a)  Martensite-austenite  phase  transitions  in  SMA.  (b)  The  residual  strain  er  for  T  <  As 
(after  [9]). 


The  shape  memory  effects  have  been  utilized  in  devices  such  as  pipe  couplers  in  which  SMA  sleeves 
are  heated  to  seal  joints  and  gas  cutoff  valves  which  close  in  response  to  elevated  temperatures.  SMA 
wires  have  also  been  employed  as  actuators  in  applications  ranging  from  remote  antenna  deployment  to 
modifying  the  camber  of  a  helicopter  blade.  Their  use  in  a  prototypical  beam  and  a  reconfigurable  wing 
are  depicted  in  Figure  28.  Shape  memory  alloys  currently  play  an  important  role  as  actuators  in  smart 
structures  requiring  the  generation  of  giant  strains  and  their  contributions  should  continue  to  increase 
as  advances  in  material  science  lead  to  improved  bandwidths. 

4.1  Material  and  System  Models 

The  modeling  of  the  constitutive  behavior  of  shape  memory  alloys  is  difficult  due  to  the  large  strains 
and  nonlinear  phase  mechanisms  inherent  to  the  materials.  The  magnitude  of  the  strains  precludes 
the  use  of  infinitesimal  elasticity  theory  and  one  must  consider  nonlinear  deformation  theory  based 
on  finite  stresses  and  strains.  As  illustrated  in  Figure  27,  the  constitutive  relations  are  also  highly 
nonlinear  and  exhibit  significant  hysteresis  due  to  the  nature  of  the  phase  transitions.  While  aspects 
of  this  behavior  have  been  quantified  through  thermodynamic  constitutive  models  [53,  56],  the  state  of 
energy-based  models  quantifying  the  dynamics  of  shape  memory  alloys  lags  behind  that  of  piezoelectric, 
electrostrictive  and  magnetostrictive  materials. 


Embedded  SMA  Layers 


(a) 


Figure  28:  (a)  Prototypical  beam  with  embedded  SMA  layers,  (b)  Embedded  SMA  actuators  in  a 
reconfigurable  wing  (after  [87]). 
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The  complexity  of  the  constitutive  behavior  and  difficulties  associated  with  inverting  energy-based 
models  to  facilitate  linear  control  design  has  prompted  the  development  of  several  phenomenological 
models  for  smart  structures  which  utilize  shape  memory  alloys.  To  illustrate,  we  summarize  the  devel¬ 
opment  of  a  Preisach  model  for  shape  memory  alloys  and  outline  the  manner  through  which  it  can  be 
employed  in  control  design. 

The  construction  of  a  Preisach  model  quantifying  the  stresses  generated  by  an  SMA  wire  in  response 
to  a  change  in  temperature  parallels  exactly  the  Preisach  formulation  detailed  in  Section  2.3  for  hysteresis 
in  piezoceramic  patches,  thus  attesting  to  the  generality  of  the  method.  In  this  case,  the  relation  between 
input  temperatures  and  generated  stresses  is  quantified  by  the  Preisach  operator  Pfl  specified  in  (15) 
where  the  input  v  —  T  is  the  temperature  and  £(s)  denotes  the  initial  state  of  the  system.  The 
Krasnolselskii-Pokrovskii  (KP)  kernel  (16)  is  employed  to  guarantee  continuity  with  respect  to  time  and 
parameters.  As  detailed  in  Section  2.3,  the  approximation  of  the  measure  p  by  a  linear  combination  of 
Dirac  measures  and  restriction  of  the  Preisach  plane  to  a  bounded  region  AS  yields  the  approximating 
operators  Pm  specified  by  (18).  The  operators  Pm  are  amenable  to  numerical  implementation  and 
inversion  to  facilitate  linear  control  design.  Further  details  regarding  the  formulation  of  Preisach  models 
for  shape  memory  alloys  are  provided  in  [6,  7,  51]. 

To  illustrate  the  use  of  the  Preisach  model  for  characterizing  the  shape  changes  produced  in  a 
structure  by  SMA  actuators,  we  consider  the  beam  depicted  in  Figure  28  which  has  a  pair  of  embedded 
SMA  layers  with  cross  sectional  area  Asm  A-  We  note  that  this  configuration  is  treated  as  a  prototypical 
example  in  [51]  and  the  following  model  is  from  that  source.  For  the  coordinate  system  specified  in 
Figure  13,  the  dynamics  of  the  beam  are  modeled  by  the  coupled  system 


d2w  d2Mint 

p-W{t'x)  +  -0£T- 

dT 


(t,x) 


d2Mext 

dx2 


(t,x,T) 


(t,  x,z)  —  k 


(  d 2 


+ 


a 2 


-  ydy2  dy2 


(58) 


T(t,y,z)  =  kv2(t ) 


with  appropriate  boundary  and  initial  conditions.  Here  w  and  T  respectively  denote  the  transverse 
beam  displacement  and  temperature.  As  noted  in  (22),  the  internal  moment  is  given  by  Mint  — 
The  density  p,  stiffness  El,  and  Kelvin- Voigt  damping  cjjl  coefficients  are  treated 
as  effective  parameters  which  incorporate  the  material  properties  of  the  beam  and  SMA  layer.  In 
the  heat  equation,  k  denotes  the  thermal  conductivity  and  v(t)  is  the  current  applied  to  the  shape 
memory  alloy.  The  external  moments  generated  by  the  shape  memory  layer  are  modeled  by  the  Preisach 
characterization 

Mext(t,x,T)  =r{x)  [  [ks{M(T),£(s)]{t)  dp(s )  (59) 

Js 

where  r(x)  is  a  shape  function  specifying  the  spatial  stress  distribution  and 

[M(T)](f)  =  — - —  f  T(t,y,z)  dydz 
Asma  Jasma 

quantifies  the  effective  temperature  in  the  layer.  Further  details  concerning  this  model  and  the  identi¬ 
fication  of  the  measure  p(s)  are  provided  in  [51]. 

For  control  implementation,  the  infinite  dimensional  models  (58)  and  (59)  must  be  approximated. 
This  can  be  accomplished  by  formulating  the  system  in  weak  form  and  employing  Galerkin  represen¬ 
tations  analogous  to  those  illustrated  in  Sections  2.4  and  3.4.  The  final  matrix  system  will  be  in  the 
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form 


y(t)  =  Ay(t)  +  [£(«)](*)  +g{t) 
2/(0)  =  2/o 


(60) 


with  the  structure  of  the  nonlinear  input  term  dictated  by  the  Preisach  model  (59). 

4.2  Adaptive  Control 

The  characterization  of  hysteresis  inherent  to  shape  memory  alloys  using  the  Preisach  model  (18)  re¬ 
quires  the  determination  of  the  weights  9  =  \ani  ■  ■  ■  ■  ,  a„m]  through  either  on-line  or  off-line  means. 
Off-line  identification  can  reduce  the  computational  overhead  for  control  implementation  but  does  not 
provide  a  mechanism  for  tuning  parameters  to  accommodate  changing  system  conditions.  This  ca¬ 
pability  can  be  provided  by  adaptive  identification  and  control  methods.  The  utilization  of  adaptive 
controllers  which  update  parameters  in  accordance  with  system  outputs  is  also  compatible  with  the 
goal  in  certain  smart  structure  applications  of  modifying  structural  properties  to  accommodate  dy¬ 
namic  exogenous  forces  or  changing  environmental  conditions.  We  outline  here  basic  mechanisms  for 
constructing  a  model  reference  adaptive  control  (MRAC)  in  the  context  of  shape  memory  alloys.  The 
reader  is  referred  to  [40]  for  details  concerning  the  underlying  theory  and  [85]  for  a  complete  discussion 
concerning  the  implementation  of  these  methods  in  a  variety  of  systems  having  actuator  and  sensor 
nonlinearities. 

The  basic  configuration  of  a  model  reference  adaptive  controller  is  illustrated  in  Figure  29.  For  fixed 
inverse  compensation,  the  inverse  B~ 1  or  approximate  inverse  B~ 1  is  determined  off-line  and  used  to 
cancel  the  hysteresis  and  nonlinearities  associated  with  B.  For  the  approximate  inverse  S_1,  the  signal 
u  to  the  plant  is  approximately  the  desired  signal  ud  specified  by  the  control  law. 

As  detailed  in  [40,  85],  the  control  ud  is  formulated  as 

ud{t)  =  9f  un(t)  +  6»2Tw2(f)  +  9ly{t)  +  0lr{t) 

where  9%  are  vector- valued  and  //[[,  9\  are  real- valued  parameters.  The  control  vectors  uq(i)  and  cu2{t) 
are  specified  by 

uui{t)  =  Ki[D)ud[t)  ,  ui2{t)  =  K2{D)y{t) 

where  K\[D)  and  K2[D )  are  filters  compatible  with  the  form  of  the  nonlinearity  (e.g.,  see  [85,  page  57]). 
Following  the  convention  employed  in  [85],  D  represents  the  time  differentiation  operator  or  Laplace 


Figure  29:  Model  reference  adaptive  control  design  for  a  linear  plant  with  a  nonlinear  actuator  and 
inverse  compensator  (after  [85]). 
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transform  variable  in  continuous  time  or  the  time  advance  operator  or  z-transform  in  discrete  time.  The 
gains  9*,  92,91, 9\  are  specified  in  a  manner  which  guarantees  that  closed  loop  signals  are  bounded  and 
tracking  errors  converge  exponentially  to  zero. 

While  fixed  inverse  compensation  is  effective  in  certain  regimes,  it  does  not  provide  the  capability 
for  tuning  the  approximate  inverse  1  to  accommodate  changing  operating  conditions  or  to  optimize 
performance.  This  motivates  the  consideration  of  adaptive  control  laws  which  systematically  update  the 
parameters  6  =  [ani ,  ■  ■  ■  ,anj J  specifying  the  inverse.  As  detailed  in  [40,  85],  this  can  be  accomplished 
through  a  variety  of  methods  including  gradient  and  Lyapunov-type  updating  laws. 

The  generality  of  this  approach  makes  it  feasible  for  several  hysteresis  models  and  their  inverses.  It 
is  illustrated  for  a  general  Preisach  model  in  [85]  and  numerical  examples  illustrating  model  reference 
control  of  SMA  modeled  using  fixed  inverses  based  on  the  KP  kernel  are  provided  in  [51].  Finally,  the 
adaptive  updating  of  parameters  in  the  KP  kernel  is  experimentally  validated  for  both  periodic  and 
random  stress  variations  in  [87].  This  latter  reference  also  illustrates  that  PI  controllers  which  ignore 
hysteresis  are  effective  at  very  low  frequencies  (<2  Hz)  but  produce  large  overshoots  and  inaccurate 
tracking  at  higher  frequencies.  This  reiterates  the  necessity  of  accommodating  the  hysteresis  in  the 
controller  to  fully  utilize  shape  memory  alloys  in  smart  structure  design. 

5  Concluding  Remarks 

The  emphasis  in  this  review  on  piezoelectric,  electrostrictive,  magnetostrictive  and  shape  memory  alloys 
reflects  their  prominence  in  current  smart  structures.  This  list  is  by  no  means  complete,  however,  and 
several  alternative  compounds  with  smart  material  capabilities  play  fundamental  roles  in  the  design  of 
structures  which  can  adapt  or  respond  to  external  stimuli.  Among  currently  employed  materials,  this 
includes  polyvinylidene  fluoride  (PVDF)  films,  electrorheological  and  magnetorheological  fluids,  and 
various  composites.  The  piezoelectric  PVDF  films  are  lightweight  and  compliant  compared  to  PZT 
patches  which  limits  their  utility  as  structural  actuators  but  provides  them  with  unique  sensor  capa¬ 
bilities  since  they  can  be  configured  to  target  specific  modes  or  structural  responses.  Details  regarding 
the  use  of  PVDF  sensors  in  structural  acoustic  systems  can  be  found  in  [17,  18].  Electrorheological 
and  magnetorheological  fluids  are  characterized  by  reversible  changes  in  the  rheology  of  the  fluids  when 
subjected  to  electric  or  magnetic  fields.  This  provides  the  capability  for  converting  the  fluids  from  low 
viscosity  states  to  ones  which  are  essentially  solid.  Potential  applications  which  utilize  these  materials 
include  vibration  isolation  systems  such  as  active  shock  absorbers  and  engine  mounts  as  well  as  electri¬ 
cally  or  magnetically  triggered  valves.  Further  applications  employing  these  materials  are  described  in 
[9,  33,  47]  while  details  concerning  the  modeling  of  electrorheological  and  magnetorheological  fluids  can 
be  respectively  found  in  [48,  66]  and  [57,  73,  74], 

As  indicated  throughout  the  review,  the  state  of  smart  structures  is  rapidly  evolving  as  new  materials, 
models  and  control  algorithms  are  developed.  To  optimize  the  capabilities  of  the  materials  in  smart 
structure  applications,  it  is  necessary  to  consider  material  characterization,  model  development  and 
control  design  in  unison.  This  permits  the  development  of  models  which  quantify  physical  mechanisms 
pertinent  to  control  implementation  and  the  design  of  control  techniques  which  capitalize  on  strengths  of 
the  smart  material  components  and  are  commensurate  with  the  models.  The  synergy  between  material 
characterization,  model  development  and  the  design  of  smart  structures  which  utilize  these  materials 
also  extends  to  the  design  of  new  smart  material  compounds.  For  example,  analysis  of  the  fundamental 
mechanisms  which  produce  the  magnetomechanical  coupling  in  Terfenol  has  led  to  the  development  of 
holmium  alloys  which  exhibit  reduced  hysteresis  [92], 

To  date,  the  majority  of  controllers  for  smart  structures  have  been  based  on  linear  models  and 
control  laws.  This  can  be  attributed  to  the  prevalence  of  PZT  transducers,  which  exhibit  quasilinear 
behavior  at  low  to  moderate  drive  levels,  as  well  as  the  relative  ease  with  which  linear  control  laws 
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can  be  implemented  as  compared  with  nonlinear  methods.  The  introduction  of  new  high  performance 
materials  and  more  stringent  demands  on  existing  materials,  however,  has  led  to  the  consideration  of 
transducers  which  are  inherently  nonlinear.  This  necessitates  the  development  of  models  and  control 
laws  which  accommodate  inherent  hysteresis  and  constitutive  nonlinearities.  The  consideration  of  non¬ 
linear  models  and  control  methods  is  also  required  for  applications  which  employ  electrorheological  and 
magnetorheological  fluids  due  to  the  highly  nonlinear  dynamics  associated  with  the  materials.  When 
characterizing  the  materials,  it  is  advantageous  to  develop  models  which  are  applicable  to  a  wide  range 
of  materials  and  can  be  inverted  (at  least  approximately)  to  construct  inverse  compensators  for  linear 
control  design.  While  nonlinear  control  laws  can  provide  optimal  performance,  the  ease  of  implementing 
linear  controllers  makes  them  preferable  in  the  majority  of  current  smart  structure  applications.  Finally, 
the  models  and  control  methods  will  continue  to  evolve  to  accommodate  the  development  of  new  smart 
materials  and  design  of  smart  structures  which  utilize  these  materials. 
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